Non-Hermitian Luttinger liquids and flux line pinning in planar superconductors 



O 

o 

(N 



O 

(N 



O 

o 



o 



Ian Affleck, ^'^ Walter Hofstetter,'^ David R. Nelson* and Ulrich SchoUwock^ 
^Department of Physics, Boston University, Boston, MA 02215, USA 
^Department of Physics & Astronomy, University of British Columbia, Vancouver, B.C., Canada, V6T IZl 
'^Physics Department, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
"^Lyman Laboratory, Harvard University, Cambridge, MA 02138, USA 
^Institute of Theoretical Physics C, RWTH Aachen, D-52056 Aachen, Germany 

(Dated: February 2, 2008) 



As a model of thermally excited flux liquids connected by a weak link, we study the effect of a 
• single line defect on vortex fllaments oriented parallel to the surface of a thin planar superconductor. 

C_J ' This problem can be mapped onto the physics of a Luttinger liquid of interacting bosons in 1 spatial 

dimension with a point impurity. When the applied magnetic field is tilted relative to the line defect, 
the corresponding quantum boson Hamiltonian is non-Hermitian. We analyze this problem using a 
combination of analytic and numerical (density matrix renormalization group) methods, uncovering 
a delicate interplay between enhancement of pinning due to Luttinger liquid effects and depinning 
P5 . due to the tilted magnetic field. Interactions allow a single columnar defect to be very effective in 

O ' suppressing vortex tilt when the Luttinger liquid parameter g <1. 

u 
u 

CLc I. INTRODUCTION 

The past fifteen years have seen much work on the statistical mechanics and dynamics of thermally excited vor- 
d tices in Type II high temperature superconductorsiSi^ii The competition between interactions, pinning and thermal 
fluctuations gives rise to a wide range of novel phenomena, including a first order melting transition of the Abrikosov 
I flux lattice into an entangled liquid of vortex filaments,^ a proposal for a highly disordered vortex glass dominated 
' O by point pinning,^ a theory of a Base glass phase with vortices strongly pinned to columnar defects^ and a distinct 
Bragg glass where point disorder converts Bragg peaks associated with crystalline order into power law singularities 

Much progress can be made on these problems using a classical continuum elastic theory^^^ An alternative approach, 
useful for treating parallel columnar defects, is to regard each vortex as an imaginary-time world line of a boson in a 
Feynman path integral, corresponding to a quantum theory of interacting bosonsj^ Here the imaginary time direction 
' T is parallel to the columnar defects. A hydrodynamic treatment of this quantum model leads naturally to the same 
, continuum elastic theory, where the classical free energy is now a classical action. If the direction of the external 
J — ' magnetic field does not coincide with that of columnar defects, it is convenient to separate the transverse component 
■r^lj- \ of the field H± from the parallel component ff|| along f. When H± << the transverse component H± plays the 
00 ■ role of a constant imaginary vector potential for the bosons. ^'^S The corresponding fictitious quantum Hamiltonian 
Op ' is non-Hermitian, with new and interesting properties. Stimulated by vortex physics, there has been considerable 
work on non-Hermitian models of non- interacting bosons in a constant imaginary vector potential, h oc H±, and a 
disordered site-diagonal pinning^iiiSii^ Since the Hamiltonian is non-Hermitian, the energy eigenvalues can be either 
■ real or complex. As discussed in |0, all states with complex eigenvalues are extended, whereas those with real 
' eigenvalues are usually localized. 
^ , Less is known about non-Hermitian models with interactionsiifi A disordered array of parallel columnar defects 
leads to a strongly pinned low temperature Bose glass phase. For h less than a critical value he, this phase exhibits a 
^ [ "transverse Meissner effect" , such that the vortex filaments remain pinned to the columns even though the external field 
O ■ is tilted away from the column direction. Although the transverse Meissner effect has now been observed in many high 
Tc superconductors with correlated disorder, recent measurementsi^ of the flux flow resistivity for h± }^ he disagree 
with a simple theory which assume that vortices tilt from column to column via kink excitations.^'"' Interactions 
' combined with vortex pinning for a periodic array of columnar defects were studied using a non-Hermitian boson 
, Hubbard model in Here, tilting the field drives a transition out of a low temperature "Mott insulator" phase 

(periodic array of vortices attached to the columns with an energy gap) into a "superfluid" phase, i.e., an entangled 
flux liquid. The corresponding non-Hermitian boson Hubbard model with site-diagonal disorder was studied in (1-1-1)- 
dimensions using a Hartree-Bogoliubov approximation in 17]. 

In this paper, we study the effect of a .single columnar pin on the statistical mechanics of thermally fluctuating 
vortex lines confined in a thin, superconducting slab. (See Fig. The external field can tilt away from the direction 

of the column, leading to statistical physics controlled by a non-Hermitian quantum Hamiltonian. As discussed below, 
the physics is equivalent to a Luttinger liquid of interacting bosons with a point impurity. Tilting the field introduces 
a constant imaginary vector potential. 

This problem is interesting for a number of reasons. If the average spacing between columnar defects is d, this is 
the regime H >> B^ = (j)o/d'^ where i?^ is the "matching field" and (f>Q w 2.07 x 10~^Gcm^ is the flux quantum. As 
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FIG. 1: Vortices with a single columnar pin (heavy vertical line) and parallel magnetic field with displacement field u{x,t) 
defined. 



emphasized by Radzihovsky,^® the large number of "interstitial" vortices between columnar defects could be locally 
crystalline or melted into a flux liquid. A third possibility is vortices in an immediate "supersolid" phase, which is 
both crystalline but nevertheless entangled, due to a finite concentration of line-like vacancy or interstitial defectSi^S 
By studying the response to a single columnar pin we can better understand the response in these phases to a dilute 
concentration of columnar defects. A dilute concentration of twin planes, a common occurance in bulk samples of 
YBCO, provides a related example of correlated pinning. 

We examine here a similar regime in (1 + 1)- dimensions, when only a single columnar pin is present. The feasibility 
of studying vortex physics in samples which are effectively (1 + l)-dimensional was demonstrated by BoUe et. a^, 
who used micromechanical oscillators to track the entry of quantized vortex filaments near Hd in a thin sample of 
the low Tc superconductor NbSe2. The observed behavior could be described in the framework of interacting vortex 
lines in a (1 + l)-dimensional random potential representing the effects of point disorder. 

Similar experiments might be possible on thin high - Tc platelet samples, with temperatures high enough to allow 
vortex interactions to screen out the effect of the point disorder. The effect of a single columnar pin might be mimicked 
by gouging a long straight scratch or notch on one side of a thin samples^ A (1 + 1) dimensional array of flux lines 
could be created by a field H, approximately parallel to the notch, and strong enough to insert a single layer of flux 
lines into a sample with a thickness given roughly by the London penetration depth. As discussed below, (1 + 1)- 
dimensional arrays of vortex lines show algebraic decay of both translational order and the boson order parameter. 
Thus, experiments and theory on this (1 + 1) -dimensional problem might give some insight into the effect of a dilute 
concentration of columnar pins on a supersolid phase in (2 + l)-dimensions. Because of the long range correlations, a 
single columnar defect can have a large effect on the flux liquid, similar to the screening cloud surrounding a Kondo 
impurity in a metal. 

A related problem in (2 + l)-dimensions concerns the effect of a single twin plane or grain boundary on vortex 
matter. A dense array of such planes parallel to the field direction leads to a Bose glass phase at low temperatures 
and a flux liquid at high temperatures?^ A single such plane should have an interesting influence on the vortex matter 
which surrounds it. Consider the effect of a small tilt on the vortex configurations. The "motion" of the tilted vortex 
configurations across the twin plane is an imaginary time version of particle transport across a Josephson junction. 
The "transport" process is likely to be quite different, depending on whether the surrounding vortex matter is in a 
flux liquid, vortex crystal or supersolid phase. A single planar defect in (2 + l)-dimensions would also have a strong 
effect on the flux flow resistivity in response to a current parallel to the plane leading to a Lorentz force which is 
approximately perpendicular to it. Motion of the Bragg planes of the Abrikosov flux lattice across such a defect is 
reminiscent of transport in materials with charge density waves,^^ 

The single defect, (1+1) dimensional physics problem is very tractible using the continuum elastic field theory ap- 
proach. Its quantum hydrodynamic formulation corresponds to a single component "Luttinger liquid" ^ The resulting 
field theory, when the field is parallel to the defect, is identical to one studied earlier^'* in the context of interacting one- 
dimensional (spinless) fermions with a point defect and of a defect in an S=l/2 Heisenberg antiferromagnetic chain.^ 
We directly apply these results to the present situation. These results depend crucially on a dimensionless parameter, 
g which controls all critical properties of the model. For instance the density correlations (with no defects) exhibit 
power law decay with an exponent rj = 2g. For interacting fermions 5 < 1 corresponds to repulsive interactions and 
g > I attractive interactions. In the bosonic case the dependence of g on microscopic parameters is more subtle, as 
we discuss. When 5 > 1 a defect is an irrelevant perturbation at long length scales, in the renormalization group 
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(RG) sense, whereas when g < 1 it is relevant, effectively flowing to infinity in the long wavelength limit. As was 
first suggested by DeGenneS)2& in the dilute limit our vortex problem is equivalent to free fermions, 5=1 pSiSLSLS^ 
We show that a transverse magnetic field defines a characteristic length scale, oc which acts as an infrared cut 
off scale in the renormalization group fiow equations. Thus even when g < 1, a defect is ultimately an irrelevant 
perturbation at sufficiently long length scales. 

We also study a tight-binding version of the non-Hcrmitian ID quantum model using the Density Matrix Renor- 
malization Group (DMRG) method generalized to treat non-Hermitian Hamiltonians. This powerful method provides 
valuable checks on our RG arguments and more detailed information about the model. 

We focus on the vortex density (i.e. magnetic field) oscillations set up by a single columnar pin as well as on 
the transverse Meissner effect. These correspond, respectively, to generalized Friedel oscillations and an imaginary 
current in the quantum problem. While these Friedel oscillations have power-law decay for a parallel field, they decay 
exponentially when h ^ Q. The imaginary current can be expressed in terms of a "pinning number", iVp, which 
measures how many bosons are stuck in the vicinity of the defect at any given "time" . We find that Np can diverge 
as — > 0, with different critical behavior in the cases g <1 and <? > 1. 

A real sample would always contains some density of point defects in addition to one or more line defect. We 
show that any finite density of point defects alters the critical behavior associated with the line defect at sufficiently 
large length scales. Thus to observe the critical behavior caused by a single line defect it will be neccessary to have 
sufficiently clean samples. 

In the next section we briefiy review the classical continuum elastic theory of fiuctuating lines and the approach 
based on mapping into a quantum model. We also discuss the value of g, showing that it goes to 1 at low densities 
and at higher densities. We note that g need not be a monotonic function of density for vortex arrays. It is possible 
that g passes through unity at some finite density as well. In Sec. Ill we study the 3=1 case, which occurs in the 
dilute limit, by exploiting the correspondence to non-interacting fermions. In Sec. IV we study general values of g 
using renormalization group and numerical methods, extending earlier results24j2£ to the non-Hermitian case. Sec. V 
contains a discussion of point defects. Sec. VI contains our conclusions. 

Appendix A derives a result for interacting bosons in one dimension of general applicability. While it is well known 
that (7 — *■ 1 in the dilute limit, corresponding to free fermions, we study the leading correction to this behavior at 
small finite density, no • We find that the result can be conveniently expressed in terms of the even channel scattering 
length, a. This quantity is determined entirely from the 2-body scattering problem, and is straightforward to calculate 
for any particular 2-body interaction. Our new general result is: 

gKl-2ano + 0[nl). (1.1) 

The scattering length, a, can be positive or negative depending on the details of the interactions, even though they 
are always assumed to be purely repulsive. 

Appendix B discusses determination of the value of g for our tight-binding model. In Appendix C we present results 
on the correlation function of the boson creation operator which is useful in confirming the RG picture regarding the 
relevance or irrelevance of a defect. Appendix D discusses the difference in ground state energy for periodic versus 
anti-periodic boundary conditions, another useful diagnostic for relevance or irrelevance of a defect. Appendix E 
contains some estimates of the effects of point disorder. Appendix F points out the connection between our model 
and one which has recently attracted attention from the string theory community. 

A brief summary of these results appeared earlier^ 



II. CLASSICAL AND QUANTUM FORMULATIONS OF INTERACTING FLUX LINES IN (1+1) 

DIMENSIONS 

Following, e.g., [Tsl l we consider a model free energy F for N fiux lines in an extreme type II superconducting 
sample of thickness in the r direction in the presence of a single columnar defect aligned in the r direction and 
located at a: = 0: 

where Xi{T) denotes the trajectory of the i*^ vortex line, Vy^^{x) is the repulsive interaction potential between flux 
lines, which can be taken to be local in r. We don't expect this locality assumption to qualitatively change the long 
distance physics in the dilute limit. The coupling measures the (attractive) interaction between a flux line and 
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the columnar defect. The tilt modulus, in the dilute limit, no A << 1, where no is the vortex density and A is the 
penetration depth, (i.e. H > Hd), for a planar sample which is invariant under rotations in the plane is 

where n — A/^ is the ratio of the London penetration depth A and the coherence length The canonical partition 
function for a system of N lines is given by the Boltzmann integral: 

^=^niIi/l?[x,(z)]e-^[{-(^»l/^, (2.3) 

where we have set kg — 1. 

From here it is possible to pass directly to a continuum elastic formulation of this model or else to use a quantum 
description where we regard Xj^r) as the trajectory of a particle and the classical Boltzmann sum as a Feynman 
path integral. We pursue the first direction in sub-section A and the second in sub-section B. Upon using a quantum 
hydrodynamics approximation to the quantum model we obtain the same continuum field theory. This field theory 
contains some parameters which we estimate, starting from the underlying vortex model, in sub-section C. 



A. Classical continuum elastic theory 

It is very convenient to pass to a continuum approximation using a coarse-grained displacement field, u(x,t), like 
that used to describe a two-dimensional smectic liquid crystal in an external field^ A standard way of doing this is 
to write the trajectory of the j vortex as: 

Xj{t) ^ jao + Uj{T), (2.4) 

(see Fig. la) where ao — is the average vortex spacing in the x-direction and define u{x,t) at the "equilibrium" 
positions of the vortices by: 

u{jaQ,T) ^ Uj{T). (2.5) 
This identification implies that the coarse grained density is; 



n{x, t) « no 



du 
dx 



(2.6) 



While this definition of u is standard also in higher dimensions, there is another way of defining it, special to 1 
dimension which has certain advantages but is essentially equivalent. We first define a field A{x, r) by the requirement 
that the position of the vortex along a constant r slice is the point x where A(x, r) = j. Thus A(x, t) can be taken 
to be a smooth monotonic function of x varying from to TV where N is the number of vortices. [An arbitrary smooth 
interpolation of A can be chosen between the points where it is integer-valued.] It thus follows that the density is: 

N „ oo 

n(x,r)^^<5[a;-x,(T)]^— ^ 0[A{x,r) - j], (2.7) 
j=i i=-oo 

where 9{x) is the step function. {9{x) = for x < and 9{x) = 1 for a; > 0.) Note that this definition implies that 
the number of vortex lines, along a constant r slice, between xi and X2 is: 

n{x,T)dx = [A(x2,t)] - [A{xi,T% (2.8) 

1 

where [A\ denotes the integer part. On the other hand, Eqs. (|2.4|l and (|2.6|l imply that this quantity is: 

n{x, T)dx « no{(a;2 — a;i) — [u{xi,t) — u{x2, t)]}. (2.9) 
Thus we write: 



A{x, t) = no[x — u{x, r)], 



(2.10) 
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and see that these two definitions of u{x,t) are equivalent in a coarse-grained limit. 
Upon using the Poisson summation formula, the density in Eq. (|2.7|l can be written: 



<^^r) = — ^ 5[A{x,t)-j]^— ^ 



dx 

] = -co 

When reexpressed in terms of u this identity becomes: 



m— — oo 



n[x, t) = HQ 

Here the reciprocal lattice vectors, {Gm} are 
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dx 



E 



27r7n/aoi m = 0,=bl. 



(2.11) 



(2.12) 



(2.13) 



We neglect the columnar pin for the moment and use the standard continuum elastic energy for a set of tilted 
vortex lines in (l-l-l)-dimensions, nameljiSiSii^ 



Fn ~ dxdr 



ic44 {dru f + ^cii {d^u f - {no(j)oH±/4:Tr) [dru) 

where cn and C44 are the compressional and tilt moduli respectively. We see from Eqs. H2.1(l and H2.4(l that: 

C44 « noei. 

Although the T = value of cn is determined by vortex interactions as 



cii -n^ j dyVi^^iy), 

thermal fluctuations and confinement entropy of interacting vortex lines lead to a value^^ 

cii cx T'^nl/ii, 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



in the limit Uq 0, where the dimensionless constant of proportionality is independent of the vortex interactions. The 
shear elastic constant cge, necessary to describe a triangular Abrikosov flux lattice in (2+1) -dimensions, is absent. 
Elastic moduli such as cn and C44 can be nonlocal (i.e., wavevector dependent in Fourier-space^), but for the large 
distance physics of interest to us here, we can take them to be constants, equal to their values on scales much larger 
than the particle spacing or the range of the interaction. 
Including the columnar pin simply adds a term to F: 

F^F-ea J drn(0,r). (2.18) 

We can include this in our elastic free energy by using Eq. H2.12|l to express n(0, r). Many long wavelength properties 
of the vortex-pin system can be studied using this continuum free energy. However, for some purposes it is convenient 
to use the quantum mechanical formulation of the model, outlined in the next sub-section. 

We note that the elastic continuum free energy of Eq. H2.14|l . with its adjustable constants C44 and en, modified as 
appropriate to account for point and/or columnar defects, is expected to give a valid description of the long distance 
physics of the vortex lattice for essentially arbitrary vortex densities. On the other hand, the simple free energy of 
Eq. H2.1|) is only valid at low vortex densities, no < A~^, corresponding to fields not to far above Hd- At higher fields 
nonlocal couplings are required to capture the physics at all length scales. Consequently, the estimate of C44 in Eq. 
(|2.15|) is only expected to be valid at low densities. 



B. Quantum formulation 



Let us return to the original discrete formulation of our free energy, Eq. H2.1|) . We may regard Xj^r), j — 1,2,3, ... N 
as the trajectories of N particles. To describe a physical sample containing vortices, the Boltzmann sum in Eq. ]2.'6\i 
could be done by first holding the entry and exit points of the N vortices at the top and bottom of the sample 
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TABLE I: Correspondence between quantities in classical vortex line problem and quantum boson problem, /i is the boson 
chemical potential, h is Planck's constant but h is the transverse field. Although the dimensions of corresponding quantites do 
not match, they do upon forming physically relevant combinations such as Ti^ /2m ^ T^/2ei. 
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fixed. The Boltzmann sum includes summing over all permutations of which vortex i, enters and exits at each of 
the prescribed entry and exit points. The 1/Nl factor is then neccesary to avoid over-counting. This expression can 
readily be seen to be the Feynman path integral for a density matrix of a system of N bosons.^ Eq. (|2.3|l can be 
rewritten in terms of the imaginary-time evolution operator g-^^^/^ as 

Z = Uf\e~^-"/^\^A , (2.19) 



where the bra and ket vectors are the initial and final states, respectively, obtained by summing over all entry and 
exit points. The quantum Hamiltonian H describes an ensemble of interacting bosons and is given by 



^^'2^l'^''^\d^~V ^+2jdxdyh{x)V{\x-y\)n{y)-e^h{Q). (2.20) 
where ^{x) is the bosonic annihilation operator and 

n{x) = i;\x)^{x), (2.21) 

is the boson number density. To account for tilting of the external field away from the direction of the pin, an 
imaginary vector potential ih has been included, thus making the Hamiltonian non-Hermitian. In the following we 
will set h — 1. We measure time, t, as well as position, x, in units of length. Thus various parameters in the quantum 
model have unusual dimensions. In Table (P) we show the correspondences between various physical quantities in the 
classical vortex and quantum boson models. Here uq = a^^ is the average number of bosons per unit length in the 
one-dimensional (ID) quantum system and w is the thickness of the slab. 

For numerical simulations it is convenient to define a lattice regularization of the model H2.20|l : 



L-lr ^ 



i=0 



H =^ -t (blb^+ie + foj+i^ie'') + —ni{n, - 1) + Vn,h,+i ~ EqUq 



(2.22) 



although this lattice has no physical meaning. Here fii = b\hi denotes the boson density on lattice site i. The 
equivalence between (|2.2U|I and 12.22|l holds for low densities per lattice site no << 1. For numerical convenience, 
the Hilbert space is restricted so that there can only be 0, 1 or 2 bosons on each site. We normally impose periodic 
boundary conditions: 

bL = bo. (2.23) 

In the low density limit, small h limit, Eq. H2.22|l reduces to the continuum model (|2.20|) with: 

t = l/2m. (2.24) 

To obtain a quantitative understanding of the analytical approximations used in this paper, and to determine 
numerically the Luttinger liquid parameter 5, we have applied a non-hermitian generalization of the density ma- 
trix renormalization group (DMRG) algorithm to the tight-binding Hamiltonian (|2.22|) . The DMRG method'^^ was 
originally invented to determine in a quasi-exact fashion numerically the properties (order parameters, correlations, 
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structure functions) of the ground or low lying excited states of one-dimensional, strongly correlated quantum Hamil- 
tonians, preferably with short-ranged interactions. In contrast to many other techniques, DMRG performance is 
typically enhanced by strong interactions. System sizes that can be studied for Heisenberg or Hubbard-type mod- 
els reach up to the order of a thousand sites. A good introduction and overview of many of the original DMRG 
applications may be found in Rcf. (S^). 

The key idea of DMRG is to grow iteratively, starting from some very small system size which can still be diagonal- 
ized exactly, a sequence of systems of linearly increasing size while carrying out reduced basis transformations at each 
growth step to keep the size of the underlying Hilbert space fixed. This reduced basis transformation is chosen such 
as to introduce the minimal error in the representation of the physical state of interest, most often the ground state. 
This is achieved by determining this state by some large sparse matrix diagonalization algorithm, and partitioning 
the entire system into two blocks A, B for which density matrices are derived by tracing out the states of the other 
block in the pure state projector, pa{b) = T^^b{A)\''P){'4'\- The eigenvalue spectra of the density matrices determine 
the new reduced bases for the system parts by choosing as new bases for the blocks a fixed number m of eigenstates 
of the density matrices characterized by the largest eigenvalues. 

This fundamental idea of the DMRG was generalized to renormalize not just quantum Hamiltonians, but also 
classical transfer matrices for statistical mechanics problems in two dimensions^ and quantum transfer matrices to 
study the thermodynamic properties of one-dimensional quantum systems^, complementing the T — results of 
the original method. Carlon, Henkel and SchoUwock'^^ have generalized the DMRG method to renormalize strongly 
non-hermitian transition matrices that originate in a master equation formulation of reaction-diffusion problems. The 
steady-state behaviour of these systems can then be derived from the left and right eigenstates corresponding to the 
eigenvalue with the smallest real part which are thus the "ground state" pair of that problem. This DMRG variant 
can be directly applied to determining the ground state eigenfunction pair of a non-hermitian quantum mechanics 
problem. 

The main numerical problem in the generalization of the DMRG to non-hermitian systems is given by the observation 
that non-hermitian diagonalization of large sparse matrices is much less stable than in the hermitian case. While 
eigenvalues can be obtained with satisfactory precision, eigenstates show small numerical inaccuracies that tend to 
accumulate during DMRG runs as these eigenstates are at the basis of the reduced basis transformations carried out in 
each DMRG step. In the unsymmetric Lanczos algorithm we have used^, the origin of these inaccuracies is mainly due 
to the inevitable loss of the global biorthonormality of the sets of Lanczos ansatz states for right and left eigenstates. 
We have achieved good numerical stability by applying a selective, very time-efficient re(bi)orthogonalization procedure 
introduced by DajiSi. 

In determining the density matrices for the non-hermitian DMRG, there is a further arbitrariness: there are both 
symmetric and nonsymmetric density matrices, defined by partial traces, Psymm — (l/2)Tr(|?/;i^)('!/;/j| -I- \'4'l) 
and Pnonsymm — T^T^l^P R.) {4' l\ respectively, where (V'lI and \iPr) are the left and right eigenstates of interest for the 
total system. The reduced basis transformations derived from the two choices are not identical. There is no clear-cut 
preference in non-hermitian DMRG: In the case of the quantum transfer matrix DMRG, which also suffers from (weak) 
non-hermiticity, the unsymmetric choice was found to be superior in precision"^*, whereas in another non-hermitian 
DMRG version, the stochastic transfer matrix DMRG, only the symmetric choice yields useful information^^''*^. 
Studies we have carried out for our problem of interest indicate that both approaches are feasible for high accuracy, 
but that numerical stability concerns favor the symmetric choice. 

With these choices made, we have studied system sizes up to L = 256 sites with up to = 64 bosons, the particle 
density being at or below a quarter. Up to m = 300 states have been kept in the reduced Hilbert spaces and found 
to give effectively converged results for currents, energies, and local quantities such as particle densities. The low 
particle density, leading to a strong arbitrariness in the insertion of particles during system growth, and the presence 
of an impurity mandate the application of the so-called finite-size DMRG algorithm, which has been applied up to 
11 times to achieve converged results. 

In applying DMRG to bosonic systems, the possibly divergent number of bosons per site has to be truncated 
algorithmically to some maximum number. As we are considering superconductors in the low flux line density limit, 
we fixed the maximum number of bosons per site to be 2. This constraint is consistent with average particle densities 
of no more than 0.25; the validity of this truncation has been checked for selected parameter sets by increasing 
the maximum number of bosons per site. It should be mentioned that these findings are not at variance with the 
statement that in Luttinger liquids with g < 1 impurities correspond to relevant perturbations and scale to infinity 
in effective field theories under renormalization group (RG) flowi24*2^ In the underlying lattice model, the resulting 
perfect pinning is effected by the generation of a very long-ranged effective local pinning potential whose strength 
decays only as a power law away from the impurity, but whose scale is essentially that of the original impurity.**^ Hence, 
we do not expect a particularly strong enhancement of the local boson density at the impurity site, as confirmed by 
our numerical results; in all runs, even at the impurity site, the boson density is well below 1 for all impurity strengths 
considered in this paper. 
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We now pass to a quantum hydrodynamic formulation of this model. ^'^ We may express the boson density operator 
in terms of a quantum field, u{x), using Eq. (|2.12|l . In order to conform to more standard notation in the quantum 
literature, we replace the displacement field, m(x,t), which has dimensions of length, with a dimensionless field, 

0{x,Ty. 

0{x,t) = -nQu{x,T). (2.25) 
Thus Eq. H2.12() for the density operator becomes: 



oo 



dx J 

We write the boson creation operator in the form: 



n{x) ^ (no + ^) J2 e2'^*™["''^+^(^)]. (2.26) 



m— — OO 



dd °" 



^^{x)o, ^no + ^ e2""["""+^(")le'*("\ (2.27) 



m— — OQ 



with, of course, the Hermitian conjugate expression for the operator, 'ip{x). Note that this formula is consistent with 
n = ^/;t^ using the fact that 5'^{x) « K5{x) where A is an ultra-violet cut-off with dimensions of wave- vector. We have 
introduced another bosonic field, (f>{x,t), which represents the phase of tp^ . We generally keep only the most relevant 
terms (in a renormalization group sense) in these expressions, writing: 

d9 

hix) ~ + — — h constant X cos{27r[nox -I- ^(x)]} 
dx 

ip\x) oc constant X e"^("^). (2.28) 

The 9 and (j) fields do not commute. In fact, in order to correctly reproduce the continuum commutation relations 
between n{x) and iIj^x), namely 

[h{x),4>^{y)]=5{x-y)^^{x), (2.29) 
[de{x)/dx, 0(y)] = -i5{x - y). (2.30) 



we require: 



Thus, we can identify: 

de/dx^% (2.31) 
as the momentum operator conjugate to 4>. Upon integrating Eq. 12.3()|l we obtain: 

[cp{x),e(y)]^{~i/2)sgn{x-y), (2.32) 

where sgn(a;) is the sign function, sgn(a;) — x/\x\. Thus we see that [dcji/dx, 9{y)] — —i5(x — y) and hence the conjugate 
momentum to 9 is 

tig = d(l)/dx. (2.33) 

We may now write a long wavelength low energy approximation to the Hamiltonian of Eq. 12.2011 (ignoring, for the 
moment, the tilt field and pinning potential) in terms of these phonon and phase variables. Keeping only terms of 
quadratic order in the fields and their derivatives, the Hamiltonian density less the chemical potential times the particle 
density may be written, for convenience, in terms of two new parameters, a phonon velocity c and a dimensionless 
"Luttinger liquid parameter" g: 



c 

H^{x) - fih{x) = - 



2 



TT \ dx ~^ 



(2.34) 



The first term comes from the kinetic energy in Eq. H2.20|l implying that 

eg no 



(2.35) 
TT m 
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The second term comes from the interaction term, which leads to: 



— « / dxV{\x\). 
9 J 

Upon canonically transforming to the Lagrangian and then going to imaginary time: r = it, we have: 

'dt = 7 ' 



and the imaginary time action 



where the Lagrangian density is 



dx / dTC^{x,T), 



C^(x,t) 



27r 



1 / d(j){x,' 



c \ dr 



f d<j>{x,T) 
\ dx 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



Alternatively, we may write H in terms of the phonon field, 6. With the help of Eqs. I|2.31|l and l|2.33|l we see that 
the Hamiltonian density can also be written as: 



Hg — = - 



Another canonical transformation gives: 



and hence the action in terms of the 6'-field reads: 



d9_gc 
9i ~ V 



Se = J j drCelXjT), 



with 



Ce{x,T) 



25 



1 f de{x,T) 

c \ dr 



d9{x,T) 
dx 



(2.40) 



(2.41) 



(2.42) 



(2.43) 



Eq. (|2.43() is, of course, just the result obtained from classical continuum elastic theory, Eq. (|2.14() . with: 

g — 
c - 

V C44 



(2.44) 



It is straightforward to add the imaginary vector potential, /i, in the quantum hydrodynamic formulation of the 
model. We begin with the Hamiltonian of Eq. H2.2()(l and then rewrite the boson creation operator, '4''^ using Eq. 
The modification of the Hamiltonian density in i/i-representation is: 



ihno d(j) 



dx 2r 



-no. 



(2.45) 



We may alternatively calculate the extra term in the Lagrangian density in 9 representation, using Eq. H2.33(l and 
then canonically transformating from Hamiltonian to Lagrangian, giving: 



Cg ^ Cg + h — . 

OT 

This is, of course, consistent with Eq. (|2.14ll . using Table I and Eq. H2.25|l . 



(2.46) 
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C. The Luttinger liquid parameter g: its significance and numerical value 



The quantity g is an essential parameter which determines all critical properties of the model, with or without a 
pin or a transverse field. In this sub-section we indicate the significance of g and estimate its value for the physical 
vortex problem in the limits of dense and dilute vortices. 

From Eq. H2.43|l we see that the density field, 9, has the correlator: 

< 6'(t, x)e{0, 0) >« ln[cV2 + x^] + constant (2.47) 

With the help of Eq. (|2.28|l , we see that the density correlation function is given by: 

^ ^ m— 1 

where the Am's are constants. Note that a conventional vortex lattice cannot form in our two-dimensional system, but 
we get instead quasi-long range vortex order with density oscillations spaced by the average inter-vortex separation, 
ttQ. The Fourier transformed density-density correlation function, or structure function, S{qx,qT), has algebraic 
singularities at the reciprocal lattice wave-vectors, Gm = 2Trmno: 

S{q.,qr)^— -^^ [forfe,9.)«(G™,0)] (2.49) 

[c [qx — Lr„i) + q^\ 

We see that S{q) diverges at q^ — ±Gi, qr = whenever g < I. See Fig. 

We can use the results above to assess the effect of a single columnar defect on the vortex density within perturbation 
theory. To lowest order in the pinning strength we find 

Sn{x,T) — < n{x,T) > —HQ 

= f dT'Cix,T-T'), (2.50) 



where 

C{x - x',T - t') =< n{x, T)n{x', t') >o - < ri{x, r) >o< n{x' , t') >o . (2-51) 

Here, < >o represents an average in an ensemble where both tilt and the defect are absent. When Eq. (|2.50|) is 
rewritten in Fourier space, the linear response to the perturbation induced by the columnar pin is determined by the 
structure function. For a system with spatial extent Lr in the time direction, we have 

Sniq-x^qr) = eoS{qx, 0) L^Sq^fl (2.52) 

The most singular response is at wave- vector {qx^qr) — (±Gi,0), i.e. for m = 1 in Eq. (|2.49|l . Thus the response 
diverges for g < 1 suggesting that a single pin is a relevant perturbation for g < 1 but is irrelevant for g > 1. This 
conclusion will be modified by a transverse field or by point disorder as we will see in later sections. The leading 
perturbative result for the oscillating part of Sn{x) at long distances from the pin in real space is thus: 

(2.53) 

Similarly, from the second of Eqs. H2.28|l and Eq. H2.39|l we obtain the correlation function of the boson creation 
operator: 

+ , , , , constant 
< .pHx,r)mO) >-> [^2 + (^,)2]i/4, - (2-54) 

This correlation function can be also derived directly from the classical continuum elastic theory^ where it is pro- 
portional to exp[— y(a;, t)/T], V{x, t) being the extra free energy arising from a dislocation pair located at (x, t) and 
(0,0). Such a configuration is shown in Fig. Q. [A similar method can also be used to explore correlations of the 
boson order parameter associated with flux lines in (2+1) dimcmsions4^] We do not usually allow such topological 
defects in our Boltzmann sums over vortex configurations. Nonetheless, we shall see in Sec. IV that this correlation 
function is very useful in studying the limit of a strong pin. 
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(a) 



S(q,,q,= 0) 




h=0 



Gj G3 



(b) 



^X 



FIG. 2: (a) Reciprocal lattice vectors (Gm,0) in the (gi,, (j'T-)-plane for /i = 0. The structure function describing vortex 
correlations in reciprocal space is singular at these points . The Fourier-transformed potential of a single attractive columnar 
defect is only nonzero on the heavy line along the ga;-axis. (b) Profile of the structure function S{qx,qT) along the axis. 
The structure function diverges near (Gi,0) for g < 1. It is also singular (although it need not diverge) at the higher order 
reciprocal lattice vectors. 




FIG. 3: topological defect configuration which occurs in the XY model but not in our vortex model. The extra free energy 
associated with this topological defect pair determines the correlation function of the boson creation operator in Eq. Ij2.54|l . 



We note, in passing, that the effect of the discrete mesh used in simulations of lattice models such as Eq. 
can also be treated in linear response theory about a continuum model. Indeed, if the average vortex separation is an 
integer p times the lattice spacing of the mesh (most of our numerical calculations are carried out for p = 4), we can 
take 



Vd{x) ~ —epCOs{2Trpx/ao) 



(2.55) 



to describe the periodic mesh to leading order in perturbation theory. Eq. (|2.55() also describes the first nontrivial 
Fourier coefficient arising from a periodic array of columnar pins. The Fourier transform Voigx) on the right-hand of 
Eq. (|2.52|l is now nonzero only for {q^, q^) = [2'np/aQ, 0) = (Gp, 0), where the structure function diverges according to 



S{qx,qr) 



[{qx 



qlt 



-')p/2 



(2.56) 



with 7]p = Igp^ . For bulk perturbations like Eq. H2.55|l . it is well known that the corresponding renormalization group 
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recursion relation readsi^ 

dep{l) 



(2-,7p/2)ep(0. (2.57) 



dl 

Hence, we conclude that the mesh is irrelevant relative to a continuum model whenever 

g>gc = 2//. (2.58) 

For p = 4, corresponding to 1/4-filling of the mesh with flux lines, gc = 1/8. Because all simulations at riQ = 1/4 in 
this paper lead to values of g significantly larger than 1/8, we can safely neglect the effect of the lattice on the large 
distance physics. 

We now turn to estimating the value of g. When the vortex lines are dense enough so that n^X ^ 1, the magnetic 
field is approximately uniform in the superconducting slab. On scales large compared to oq and A, we then expect 
that the energy can be approximated by the usual form from magnetostatics 



B^(r) H 



B(r) 



(2.59) 



where B(r) is the magnetic field intensity, H is the applied field and the integral runs over x, r and the short 
direction y of a thin slab of thickness w ^ \. We take H parallel to the time-like direction r and note that Eq. (|2.59|) 
is minimized for B(r) — Ht = Bqt. We then neglect the variation of B(r) across the slab along y and expand 
about the state of a uniform field by setting B(r) = {Bx,0,Bo + 6 Br), with \Bx\, \ SBr\ ^ Bq. Upon making the 
identifications SBt-{x,t)/ Bq = — 9^u(a;,r) and B.j.{x,t) / Bq ~ 9^7i(a;,r), we obtain an expression like Eq. H2.14|l with 

wB?{ wBl 
C44 = -^, cn = ^ (2.60) 

After setting Bq sa (^qUq/w, where (j)Q = 27r?ic/2e w 2 x lO^^gauss-cm^ is the flux quantum, we find from Eq. (|2.44|) 
that 

g~—^ (2.61) 

With w w lO^^cm (we have in mind thin slabs with thickness of order the London penetration depth), and T w 77 K 
we find that g w 10^'^ ^ 1 in this dense limit. 

In the dilute limit, uq ^ 0, g ^ 1, corresponding to free fermionsiS2iSS*2LS2i That 5 = 1 corresponds to free fermions 
can be simply checked from the fact that the density-density correlation function in Eq. (|2.48|) reduces to that of free 
fermions in this case: 



< n{x, r)n(0, 0) >^ nl + 



—r + A2 cos(27moa;) 



? ? 2 
c r — X 

{c^t'^ + x^y 



(2.62) 



Note that this asymptotic behavior of g is consistent with the behavior of C44 and cn in Eq. H2.15|) and H2.17|l . Free 
fermion behavior in the dilute limit follows from the equivalence of hard-core bosons with free fermions in (1-f 1) 
dimensions. As long as the average boson separation is large compared to the range of the inter-boson interaction, 
this "hard-core" result holds at long distances even for a finite range inter-vortex interaction, g = 1 is the border-line 
case where the pin is marginal. We will exploit this equivalence of dilute bosons to non-interacting fermions in the 
next section. 

In Appendix A we derive the leading correction to this result at low density, Eq. Hl.l(l . We expect this result to 
be exact for our vortex system despite the fact that we ignored vortex-vortex interactions that were non-local in r in 
our approximate free energy of Eq. H2.1|l . Such non-local interactions ultimately have similar effects to quartic and 
higher terms in the dispersion relation (see Appendix B of the second article in Q) and we expect that they only 
affect g at ©(tt-q). 

Whether g increases or decreases as the density is increased from depends on the sign of the scattering length, 
a, for the vortex- vortex interaction potential V{x). This sign depends on the detailed form of V{x). [Note, for 
example, that an infinite hard core repulsion gives a > whereas a repulsive (5-function potential gives a < 0.] The 
appropriate potential, V{x), is determined not only by the bulk inter- vortex interactions but also by the thin-slab 
geometry. Results on this scattering length will be reported elsewhere.*^ If a > then g may decrease monotonically 
with increasing no so that the pin will be relevant for any no- On the other hand, if a < 0, then g initially increases 
with no- Since we have shown that it goes to a very small value at large densities it must then exhibit non-monotonic 
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behavior. Furthermore, it must pass through the value g = 1 at some finite critical density, Uc. In this case, the pin 
will be irrelevant for uq < and relevant for rig > ^c- 

The value of g for our lattice model as a function of the microscopic parameters U/t, V/t and uq can be determined 
at low riQ from an exact formula for the scattering length, a, of the microscopic model or more generally by numerical 
methods. See Appendix B. 

III. DILUTE LIMIT: FREE FERMIONS 

In this section, we consider the dilute limit where we may approximate interacting bosons by non-interacting 
fermions. This regime corresponds to 5 = 1: a marginal pin. 

1. Density oscillations 

The density oscillations induced by the pin can be calculated straightforwardly using the standard method of 
calculating Friedel oscillations for non-interacting fermions, suitably generalized to the non-Hermitian case: 

<0|n(a;)|0>=5]V^(a:)V^(a;). (3.1) 

n 

Here and are the left and right cigcnfunctions of the non-Hermitian single-body Hamiltonian of Eq. (|2.20() 
with V set to zero. The sum is restricted to the set of single particle levels with Re{En) < Ep, where En are the 
single particle levels and Ep is the Fermi energy. (We denote this restricted sum by ^'.) 

The exact single particle energy eigenvalues and eigenfunctions for a 5-function potential with an imaginary vector 
potential were obtained by Hatano and Nelson (H-N)ii2i [See the Appendix of the second reference in (10).] Here we 
adapt these results to obtain the Friedel oscillations. Note the change of notation from that paper: 

h ^ 1 
g h 

Vo Eq 

E. ^ L. (3.2) 

We must also keep in mind that, assuming periodic boundary conditions on the bosons, the boundary conditions on 
the single-particle wave-functions in the effective fermion problem are periodic if the total number of bosons is odd 
but anti-periodic if the total number of bosons is even. [See Appendix A.] As H-N show, the (5-function potential 
produces a single bound state if meg > h. There are also extended states which are responsible for the long distance 
density oscillations. 

We consider the non-Hermitian single-particle Schroedinger equation for the right eigenfunctions: 

ilj^-{x)=Ei!^-{x). (3.3) 

The extended right eigenfunctions are written in terms of complex right- wave- vectors, X„: 

^^{x) = A„e(*^"+'')" + B„e(-*^"+'')% {x ^ 0) (3.4) 

for constants An and _B„ determined by the boundary conditions. Assuming N odd, so that the eigenfunctions obey 
periodic boundary conditions, the allowed values of K satisfy: 

K[cosh{Lh) - cos{LK)] + meosini/f = 0, (3.5) 

where we may assume, without loss of generality that ImK > 0. Provided h > and hL >> 1, Eq. (|3.5() can be 
approximated as: 

K [e^^ - e-'^^] + imeoe-'^'^ = 0, (3.6) 

or 

{K - imeo)e-'^^ = Ke^^. (3.7) 
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We see that ImK « h and hence: 



Kn ~ kn + ih + 



icj){K) 



where n runs over all integers and 



kn = 2Tm/L, [h > meg) 

kn = 2n{n + l/2)/L, {h < meo) 



(3.8) 



(3.9) 



(We may assume |Im0| < tt.) Substituting in Eq. 13.7|l implies: 

kn + ih 



HKn 



In 
In 



kn + i{h — meg) 
kn + ih 



The corresponding energies are: 



kn + i{h — meo) 



E = {l/2m)K^ 



, {h > meo) 

{h < meo). 



(3.10) 



(3.11) 



The exact right eigenfunctions, for periodic boundary conditions, are given in (A. 13) of H-N. Taking the large L limit 
these become: 



x>0. 



(3.12) 



Here, the + or — minus sign apply to h > meo Sind h < meo respectively. For A: > Eq. H3.12|) has the interpretation 
of a particle coming in from the left and being reflected and transmitted. Note that the reflected wave decays 
exponentially as e^''^ {x < 0) unlike in the normal, Hermitian, case. The corresponding left eigenfunction is given by: 



ijt{x;k)^tP\{x;k)*. 

It follows from Eq. I|3.5|l that Kn{—h)* = Kn(h). The left eigenfunctions, at large L, are given by: 



-ikx 



+ (±e'^-l)e^'=^-2/»a;^ 



(3.13) 



(3.14) 



For fc > 0, Eq. H3.14|l describes a particle arriving from the right and being transmitted or reflected, with the reflected 
wave decaying exponentially. The product of these eigenfunctions, at large L, (for either sign of x) is given by: 



4,^{x)i}^{x) = ±e^[l + {±e^ - i)e2(''=-'')l-l] 



(3.15) 



Upon normalizing these wave-fuctions and then integrating over all levels below the Fermi surface, we find the density 
at large |a;|A;j;^: 



n{x) 



kp 

TT 



-2h\x\ 



dk 

2^ 



(±e 



0(fe) 



l)e 



2ik\x\ 



TT 



ATTilxl 



.e^^k^lM + c.c. 



meo 



(h — iTTUo — meo)47ri|a;| 



e ' ' + C.C. 



(3.16) 



Note that the Friedel oscillations decay exponentially with a decay length given by l/2h. This result arises from the 
exponential decay of the reflected wave which leads to the exponential decay of the interference between incident and 
reflected waves. 

In the lattice model (|2.22l) the limit g = 1 corresponds to U oo, y = and thus to a system of hardcore bosons, 
which in (1+1) dimensions can be mapped onto the Hamiltonian for non-interacting fermions (described by creation 
and annihilation operators c\ and Ci), namely 



H = -eono-t^{c\ci+ie + cj+ice'') 



(3.17) 
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where no = CqCq and we assume a canonical ensemble with TV particles. This Hamiltonian can easily be diagonalized 
exactly for system sizes of a few hundred lattice sites. Using Eq. 1)3. l|l . we have obtained numerical results for the 
density profile shown in fig. 01 In agreement with the analytical predictions above we find strong Friedel oscillations 
which are suppressed drastically by a finite nonhermitian term h > 0. We have also confirmed that the oscillations 
decay as exp{—x/^±) with a decay length ~ 1/h. 

In fig. we give a schematic representation of the underlying "traffic jam" picture. Since flux lines repell each 
other, they form a local vortex lattice, aligned with the imaginary time direction, in the vicinity of the pin. Finite tilt 
destroys this effect at large length scales, with a crossover scale given by £^±. In fig. |Slwe have arbitrarily chosen to 
show a queue which is two deep on either side of this (symmetric) traffic jam, corresponding to a total of five maxima 
in a density plot such as Fig. 10}. 




0.22 - 

1 I 1 I 1 H I H , I , I , 

-30 -20 -10 10 20 30 

X 



FIG. 4; Friedel oscillations of the flux line density vs. distance x from the defect, calculated numerically in the free fermion 
limit {g — 1) with eo = 2 and no = 0.25. 




X 



FIG. 5: Schematic picture of the "traffic jam" scenario, for vortex lines in the vicinity of a columnar pin as described in the 
text. Because of their mutual repulsion, flux lines queue up in the vicinity of the pin. 
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2. Current and Pinning Number 



Conservation of "charge", i.e., the number of flux hncs in a slab, imphes the existence of a conserved current 
operator, even in the presence of a pin and a tih fleld. This is given by: 



J{x) 



dn{x 
dh 



—I 

2m 



dx 



dx 



In the absence of a pin there is an imaginary current in the ground state, 

ih 

J =< 0| J(a;)|0 >= — < 0|7/;^(a;)V'(x)|0 >= ihno/m. 



(3.18) 



(3.19) 



Physically, this current describes vortex lines tilted at an angle of tan~^ (/i/m) relative to the r-axis. Equivalently, 
each boson has an average imaginary-time "velocity" of h/m. Intuitively we might expect that the pin would reduce 
this current since vortices tend to get "stuck" on it. However, the current is conserved and further describes a r- 
independent vortex density, even in the presence of a pin. Thus it has the same value in the vicinity of the pin as 
everywhere else in the sample. Although a single pin cannot change the value of the current in the limit L — > oo, this 
defect nevertheless can have important finite size effects on the current which we now consider. We find it convenient 
to define a "pinning number", Np, which describes the finite size reduction of the current due to the pin. Setting 
J =< 0| J(a;)|0 >= {ih/mL){N - Np), thus defining 



Np EE N[J{Q) - J(eo)]/J(0) = N + imLJ{eo)/h. 



(3.20) 



Since we may think of each boson as contributing ih/mL to the current in the clean system, Np measures the effective 
number of bosons which are not contributing to the current because they are "stuck" in the vicinity of the pin. We 
are considering a finite size effect since we expect that J(eo) — J(0) cx 1/L. 

Np may be readily calculated in the dilute limit where we can use the free fermion approximation. The simplest 
procedure is to calculate the ground state energy and then differentiate to get the current using: 



J 



—i dEo 
T~dh' 



(3.21) 



which follows from Eq. H3.18|l . The ground state energy, Eo{h) may be calculated by summing up all single particle 
levels below the "Fermi surface" , as indicated by the prime in the summation: 



En = 



^ 2r 



(3.22) 



where the are given by Eqs. (|3.8|l and (|3.1U|) for the extended states. 

We will assume that the limit L ^ oo is being taken in Eq. (|3.2UI) so that we only need the energy to 0(1/L). In 
this way we obtain: 



ihrin 



2'KmL 



kp In 



kj^ + {h- meoY 



+ ih In 



(kp + ih)[kp — i{h — meo)] 
{kp — ih)[kp + i{h — men)] 



(Here kp = TmQ.) Upon extracting the pinning number, Np, from Eq. (|3.23|) . we have: 



kpln 



k^ 



kp + {h- meoy 



+ ih In 



{kp + ih)(kp — i{h — mep)) 
(kp — ih){kp + iih — mep)) 



(3.23) 



(3.24) 



Despite the fact that there is a bound state for h < meo but not for h > meo, is a smooth function of h near meg 
in this large L limit. (However, a step develops at /i = meo in the limit kp — > 0.) In the limit h << imo, we find the 
asymptotic behaviors at large and small eo^ 



N„ 



{meo << h) 



{meof 
2Tr^noh ' 

^ln(|eo|/»^o), {meo»h,no). 



(3.25) 
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Remarkably, the pinning number diverges as /i — > 0. This behavior can be understood in terms of the Friedel 
oscillations discussed in the previous sub-section. These oscillations imply a local density wave which extends out 
to a distance of 0{l/h) away from the pin as illustrated in fig. El We may think of the particles as entering a sort 
of "traffic jam" near the pin similar to the one which may occur near a toll booth. Heuristically, we think of each 
particle as waiting at the locations of the peaks in the density, which have spacing l/riQ, until the particle in front 
has moved ahead one space. The pin at a; = corresponds to the toll booth. Unlike most real traffic jams at toll 
booths, this one is symmetric under x —x. The number of particles participating in the traffic jam is 0{l/h) and 
represents the number of particles which are not participating in the current. 

It is important to note here that we have taken the limit L oo first before considering small h. That is, we are 
assuming L » 1/h. The behavior of the current as h —^ with fixed L is quite different, becoming linear in h. In 
this linear response regime the non-Hermiticity of the Hamiltonian becomes unimportant since the linear response to 
the imaginary vector potential can be expressed in terms of the susceptibility calculated ai h — 0. Apart from a factor 
of i, this is the same as the linear response to a real vector potential. A constant real vector potential of strength 
~ih corresponds to a dimensionless magnetic flux, ihL threading the ring. The resulting real current is known as 
the persistent current. It has been discussed for the case of a localized impurity potential in a Luttinger liquid, by 
Gogolin and Prokof'ev.^^ In the case of non-interacting fermions, these authors find an expression for the persistent 
current in terms of the transmission coefficient at the Fermi surface, Tp. For the particular case of the (5- function 
potential of Eq. H2.20|l . the transmission coefficient is given by: 

The Gogolin-Prokof 'ev formula for the persistent current, in the limit where the flux goes to zero, for the case of N 
odd [Eq. (1) of Ref. ^] then gives: 

J^^^^tan-ifl^V (3.27) 

Note that this result, linear in h, is independent of L. Naturally, at eo 0, it reduces to our previous result for the 
system with no impurity. Of course, once we go beyond linear order the dependence of the current on a (real) flux is 
very different than its dependence on an imaginary vector potential h. In particular the flux dependence is periodic 
with period 2tt and is 0{1/L). We expect the imaginary current to cross over from the large L result to the linear 
response regime when 1/h is of order L. At this point the "traffic jam" is filling the entire system. 

In the limit eo oo, tunnelling of particles past the pin becomes very ineffective. It is then instructive to rederive 
our results using a weak tunnelling model. This approach will be very useful in Sec. IV when we consider the case 
g < 1. In the dilute limit we can again analyse a non-interacting fermion problem. It is convenient to consider a 
non-interacting fermionic tight-binding model with a weak link between sites L and 1: 

L 

H = -<^[e-''4c,+i +6^1+1^] -r[e-''c[ci +e''cIcL], (3.28) 

2=1 

where T << t represents the weak link caused by an impurity with very large eo- Here c,; is a fermion annihilation 
operator. For some purposes it is more convenient to make a similarity transformation to a different Hamiltonian 
which has the same (right and left) eigenvalues, chosen so that all the non-Hermiticity resides on the weak link. This 
transformation is equivalent to the replacement: 

c] -> c]e-''J'. (3.29) 
This non-unitary, commutation-relation preserving transformation changes the Hamiltonian to: 

L 

H = -tY^[clc,+i + h.c] - r[e-''^c{ci + e^^c\cL\. (3.30) 

i=l 

The lattice Schroedinger equation associated with H is: 

-t(i/-j_i+^,+i) = Eij;, {j^l,L) 

- Fe'^^Vi = E^L. (3.31) 
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To find the scattering states we use tlie ansatz: 

= Ae'^^ + Be"'^^, (3.32) 

wliere A and B are amplitudes and K will, in general, be complex. Without loss of generality, we may assume 
ImK > 0. The eigenvalues are then: 

E{K) = -2t cos K. (3.33) 
The last two equations of H3.31|l can be rewritten as: 

Te-^^ {Ae''^^ + Be^'^^) = t{A + B) 

Te''^ {Ae'^ + Be-'^) = i Ue'^(^+i) + Be-'^"(^+i)j (3.34) 

Upon solving for A/B and simplifying we find an equation which determines the eigenvalues, i.e. the K^s: 

- 2rt cosh{Lh) sin - sin (i - 1) + sin + 1) = 0. (3.35) 

For large Lh, we see that K must have the form: 

Kn = ih + 2TTn/L + i<j>{Kn) /L, (3.36) 

as before. [We are interested in the small F case, where a bound state occurs, so Kn is given by the first line in Eq. 
(|3.8() .] After dropping terms suppressed by e"''^ (we assume h > 0), we find: 

- Ft sm{ih + k) + —e-^+'^+'l' (J^'^k+'P ^ 0. (3.37) 

2i 2i 

Here we have set k = 2TTn/N and dropped the subscript n. We can now determine as a function of k. 
Once we have 0(fc), we can calculate the ground state energy and hence the current using: 

£■0 = -2t^cosi^. (3.38) 

K 

Here we again sum over all energies whose real parts lie below the Fermi surface. 

Let us now just focus on the small F limit. In this limit we may drop the second term from Eq. (|3.37|) . It is 
interesting to note that this approximation corresponds to dropping the second term from Eq. (|3.35|) . which in turn, 
corresponds to setting Te~^^ to zero in the first of Eq. (|3.34|l . Thus we consider a "one way" model which ignores 
the hopping from 1 to L but allows it from L to 1. Eq. (|3.37|) then reduces to: 

Fe-^''sinif==sini^(i + l). (3.39) 

Note that, in this small F limit, F and h only appear in the combination Te^^. It is therefore convenient to define a 
shifted h variable: 

and define 0' by Eq. H3.36|l with h replaced by h' . Thus: 

0' = (/)-lnF. (3.41) 

The dominant correction to the current, for small F, is given by using the formula for the current with no pin, and 
replacing ft, by /i'. The redefined phase shift is now determined by: 

e-^' = l-e^''-^. (3.42) 

At this point it is convenient to take the contimmm limit, assuming that \h + ik\ << 1, so that the phase shifts 
become: 

0' = ln{-2iK), 

(j) = lTii-2iKT). (3.43) 
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Upon comparing to Eq. (|3.10|l . we see that Eq. H3.43|l is the same formula for the phase shift obtained from the 
original Hamiltonian in the limit of large Cq, with the replacement: l/meg — ^ 2r. The corresponding cm'rent is: 



ih'no iuQ 

J = 1 ln(nomj. 

m niL 



Upon using Eq. 



our final result for the current becomes: 



J 



mo 
mL 



ln(rnor7i). 



(3.44) 



(3.45) 



Because this analysis applies in the limit F ^ we see that Eq. (|3.45|l represents a reduction of the current. The 
leading dependence on T follows immediately from the formula for the current with no pin upon the replacement of 
h by h' defined by Eq. JTini : 



h' = h 



InT 



The corresponding pinning number is: 



Np^N + imLJiV) /h ■ 



no 
h 



InT. 



(3.46) 



(3.47) 



This diverges logarithmically, as F — > 0, i.e. as the pinning strength goes to infinity. It agrees with the large eo result, 
Eq. (|3.25|) . with the identification F = no/eQ. 

We have also calculated the pinning number via exact diagonalization of the noninteracting fermion tight -binding 
model H3.17|) . [See Fig. ©.] As expected from H3.25|l . Np in fact grows, with decreasing h, as 1/h until at h (x 1/L 
this divergence is cut off due to finite-size effects. As a result, in the linear response limit Lh — > the pinning number 
saturates at a value Np (x N where N is the total number of bosons. Another prominent feature is a steplike decrease 
of Np close to he = meo, which is a vestige of the single vortex depinning transitioui^ After being smeared out by 
interactions this step is only visible at low filling. 




FIG. 6: Pinning number in the free fermion limit (g — 1) for L 
h and the "step" at he ~ meo due to single-vortex depinning. 



100 and eo = 2. Note the strong enhancement at small tilt 



IV. RENORMALIZATION GROUP APPROACH TO GENERAL g 



As noted in Sec. IIC, when h = 0, a, pin is a relevant perturbation for g < I and irrelevant for g > I. We use the 
hydrodynamic approach, involving the dimensionless displacement field, 9{x, r), defined in Sec. JIB. Upon introducing 
the impurity via Eq. H2.18|l . keeping only the most relevant parts of ^(O), from Eq. p. 12(1 . and including the tilt field 
from Eq. ifOHjl . the La grangian density becomes: 



Co 



TT 

2^ 



1 rrn 

c Kdr 



de 

dx 



hdri 



eo-;- + eo cos(27r6^) 
ax 



5{x) 



(4.1) 
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where eg cx eg, the pinning strength, but includes effects of ehminating short distance modes. The lowest order 
renormalization group scaling equations for these 2 boundary interactions are well known. We study the scaling of 
the dimensionless quantities cq and eg, considering the effect of eliminating short distance degrees of freedom of 6*, 
thus reducing the short distance cut-off A. eg is unchanged by this reduction, indicating that it is a marginal coupling 
constant. In fact, we may eliminate this term completely by the transformation: 

eix) eix) + ^sgn{x) (4.2) 
27rc 

where sgn{x) is the sign function with sgn{Q) = 0. This transformation has no effect on the cos[27r0(O))] or dO/dr 
terms in Eq. (|4.1|) . This transformation, Eq. (|4.2|l . shifts the density of bosons by {eog/T^c)5{x). Thus, the eo term 
in Eq. H4.1() is an exactly marginal interaction. On the other hand, the eo term has a g-dependent scaling equation. 



A. h=0 case 



We can determine the scaling equation for eg by observing that if we integrate out Fourier modes of 9 with wave- 
vectors between Aq and A then cosM(a:,T) gets replaced by: 



A 



a 



cos27r6l 1^— j cos 2710, (4.3) 

implying that the operator cos 27r0 has a renormalization group scaling dimension of g. Noting that the eo term in 
the action involves a r-integral but no x-integral, due to the 5{x) factor, we see that 



Equivalently we may write the RG scaling equation: 

^eo-(l-5)eo, (4.5) 

where / = ln(Ao/A). The dimensionless pinning strength gets larger at longer length scales for g < 1, but gets smaller 
for g > 1. 

This model, Eq. (|4.1|l . has been well-studied in the closely related context of a quantum fcrmion system^'' or a 
quantum spin chain^^ with a point impurity, from which it arises by bosonization. In these contexts it has been rather 
well established by numerical and analytic work that, for g < 1, starting even with a small impurity strength the long 
distance (low energy) behavior is that of a "cut chain" with the large impurity strength effectively decoupling the two 
sides. On the other hand, for g > 1, starting even with a large impurity strength the long distance behavior is that of 
a "healed" chain with no impurity. We believe that this is also the case for the bosonic version of the model, defined 
by Eq. ^T^. 

The cut chain fixed point is easily studied in the phase boson representation. A boundary condition, 6(Q) = must 
be imposed. Note that we should really think of the cc > and x < Q parts of the systems as being independent in 
this limit. (We take i — > cx3 for this discussion.) Thus we get two boundary conditions, 0(0^) = 0. These imply 
89 /dt cx lie = and hence, from Eq. (|2.33l) . the Neumann boundary condition on </>: 

f(«-)^0. («) 

This boundary condition modifies the correlation functions. Of course any correlation function of two fields on opposite 
sides of the pin is zero. The correlation function of two fields on the same side is also modified. One way of calculating 
these correlation functions, with the boundary condition, is to decompose the free boson fields <j) into left and right 
moving components: 

(t){x,t)=(j)L{t + x/c)+(i)R[t~x/c) (4.7) 

(The most general solution of the equations of motion, [d]; — c^d1)(j) = can be decomposed in this way.) Eqs. (|2.33|l 
and H2.41|l then imply: 

9{x,t) = ^[4>L{t + x/c)-Mt-x/c)]- (4.8) 

TT 



21 



The boundary condition can thus be written: 

M-t,0^)^<t>B.{t,0^). (4.9) 

Let us focus on correlations for a; > 0, for example. Then the boundary condition of Eq. (|4.9f) implies, since 4>r^l is 
a function of t =p x/c only, that we may regard 4>r{x) as the analytic continuation of 4>l{x) to the negative axis: 

= 0L(-a;), (fora;>0). (4.10) 

This has the effect of making the density and boson creation operators bi-local: 



noe 

n{x) — no sa —^^^ + constant x cos{27rnoa; + 2(7[(/)l(x) — 0l(—x)]}. (4-11) 
TT ax 

The correlation functions can now be calculated using: 

< (j)L(x, T)(j)L(y, 0) >= - — ln[(a:; - y) + icr] + constant (4-12) 

4.9 

Thus the correlation function of the boson creation operator, discussed in Sec. IIC, becomes: 

f xy 1 ^^^^ 

< V^(x,r)V^(y,0) >oc , (4.13) 

where x and y are on the same side of the pin. Note that in the limit x,y » \x — y\, c\t\ we recover the bulk behavior 
[Eq. P31jl ]: 



<^^x,T)^{y,0)>^\- , , ,A . (4.14) 



1/45 

[{x - yY + c^T^]" J 

On the other hand, in the limit c\t\ >> x,y we obtain the "boundary critical behavior": 



1 1 ^'^ 

<i/.t(x,T)V'(2/,0) >(x <j ^ ^ . (4.15) 



Thus we see that the operator, -0^, which has a bulk scaling dimension of 1/4(7 ha-s a boundary scaling dimension 
which is twice as big, l/2g. To understand this result, note that, without the boundary condition, (p^ and (p^ are 
independent fields, so that both factors e"^^ and e*"^'* contribute equal amounts l/8g to the scaling dimension of . 
After imposing the boundary condition becomes the operator e^*'^^^°-' which has dimension l/2g. 

When the pin is relevant, g < 1, we may calculate the density oscillations at long distances from the pin by assuming 
that eo — > oo and using the Neumann boundary condition of Eq. H4.6|) . This constraint leads to Eq. (|4.11|) which 
leads to: 

, , , , constant • cos(27rnnx) 
<0|n(a;)|0>->no + (4.16) 

For an irrelevant pin, we expect the result of lowest order perturbation theory in e to be valid at long distances, 
giving Eq. H2.53|l . Note this involves a different (larger) exponent than the one which occurs for a relevant pin. These 
density oscillations are the bosonic version of the generalized Freidel oscillations discussed for fermionic systems and 
spin chains in Ref. 

We can now study the stability of the cut chain fixed point. As discussed above, in this limit the system decouples 
into two separate sections to the left and right of the pin. If eg is very large but finite, there will be a (dimensionless) 
weak tunnelling matrix element, F cx 1/eoi between the two sides. In the ■0 representation, this effective Hamiltonian 
is: 

H = + H+ - Tc[^p\0~)iP{0+) + h.c] (4.17) 

Here 

H+=[ dx-^^^ + constants ipUO)ip{0) + ^ [ dxdyn(x)V (\x ~ y\)h{y) , (4.18) 
Jx>o 2m dx dx 2 J^^y^o 
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iJ_ is defined similarly. (Alternatively, an interaction term similar to that in Eq. (|2.22() may be maintained between 
the two sides. The important thing is that there is no motion of bosons across the pin in this limit.) This is 
conveniently written in terms of the phase boson, cf) with the boundary condition d(t)/dx = at x = 0*. Thus we get 
two copies of the Lagrangian of Eq. H2.39|l for a; > and x < 0. As explained above, the tunnelling term, analogous 
to a Josephson coupling across the weak link, becomes: 

- Tct/;^ (0^)^(0+) + h.c. cx rcos2[0L(O+) - (/'l(0")], (4.19) 

of scaling dimension l/g. (Because "0(0^) s-iid V'(0 ) become independent operators upon imposing the boundary 
condition, their scaling dimensions simply add. ) Thus, the renormalization group equation obeyed by F is: 

We conclude that a weak tunnelling across the pin is relevant for 5 > 1 but irrelevant for g < 1, i.e. precisely the 
inverse of the situation for a weak pinning potential, eq. This implies consistency of the bold assumptions that a weak 
eo will renormalize all the way to the broken chain fixed point (corresponding to F ^ 0) for g < 1 and that even a 
large eo will renormalize to for g > 1. We present numerical evidence to verify this conjecture in Appendix C and 
D. 



B. Density oscillations for h > 

We now consider the effect of a non-zero transverse field, h, in Eq. H4.1() . We may eliminate the term in Eq. 14. 1() 
proportional to ft. by a shift of the 9 field: 

0{x, t) 0{x, t) ~ igch/TT)T. (4.21) 

Let us first consider the density oscillations in the limit of a weak pin, in lowest order perturbation theory in eo as in 
Sec. lie. The singular part of the density-density correlation function picks up an extra /i-dependent phase from this 
shift: 

, , qc^ c^t"^ — x'^ coa\{2Tmnx — 2qhcT)p\ . . „„x 

< n{x, r n >^ + , + \\ \ J ^ + . . . , 4.22 

Upon Fourier transforming, we see that the singularities of the structure function have moved off the qx axis to 
(Qk, It) — (27mo, 2ghc)p, as shown in Fig. (|7|). The linear response to a static pin is proportional to S{qxj 0), [see Eq. 
(|2.52(l ]. which is non-singular, as indicated in Fig. Focussing on what was the leading singularity, at q^ = ^tthq, 
we find: 

dr-^ . (4.23) 



This integral may be expressed in terms of the modified Bessel function of the second kind, A'g_i/2 

^i(2iThno / m)T 2 / 

■ = ( — 

[3.2 _,_ ^2^2]g ^^2g^lY{^g) \2\X 



I 



„i{2TThna/m)T 9 / f \ 9^1/2 

^^TzmriT^- 772^^17^ Un:T ^^-1/2(^/0, (4.24) 



with characteristic length scale 

C = cm/(27rno/i) = l/(2.g/i) (4.25) 

and where T{g) is Euler's F function. In the limit \x\ << wc obtain the zero tilt result, Eq. 12.531) . In the opposite 
limit, |a;| >> ^, we obtain exponentially screened Friedel oscillations: 

< 0|n(a;)|0 >« no + constant • eo • cos(27moa;)e"l^l/^-— p. (4.26) 

Beyond the new characteristic length scale, ^ = l/{2gh), we expect that the pin loses its effectiveness in ordering 
the vortices, for any value of g. It is natural to assume that this new length scale acts as an infrared cut-off on the 
renormalization of the pinning strength, e{l). Thus the long distance physics should be controlled by: 

eo{h) « ~eoo{Ao/h)^~s. (4.27) 
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(a) (b) 

FIG. 7: (a) Reciprocal lattice vectors in the (g^;, (j'T-)-plane for h ^ 0. In contrast to Fig. J^J, the singularities now occur off 
the qa;-axis, where the potential due to the columnar pin is zero, (b) As shown in this profile, the structure function S{qx,0) 
is now finite and nonsingular everywhere along the q^-axis. 



Here eoo is the "bare" value of eg, that is, before renormahzation. Provided that eo(^) << 1, we expect the perturbative 
result to be valid at long distances. Thus it would be valid for g > 1 even if the bare dimensionless pinning potential 
is not small, for sufficiently small h. It would also be valid for g < 1 provided that the bare dimensionless pinning 
potential is sufficiently small and h is not too small. As we showed in the previous section, for the case g — 1, 
corresponding to free fermions, this exponential decay of the Friedel oscillations at long distances holds for any eo- 
Even if the above criteria are not satisfied we still expect exponentially decaying density oscillations at long distances 
for any values of g and eg provided that h ^ Q. 

DMRG results for the Friedel oscillations with g < \ are shown in fig.|Hl Note that, as in the free Fermion limit g — 1, 
a finite tilt h tends to strongly suppress the oscillations. Our data are consistent with a crossover from power-law to 
exponential decay with increasing /i, although the system sizes are too small to extract precise exponents. 
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FIG. 8: Friedel oscillations calculated by DMRG for eo = 2, C/ = 10, F = 4, L = 128 and no = 0.25, corresponding to 5 « 0.72. 




C. Current 



We first consider the case h ^ for fixed L. As discussed in Sec. Ill, when hL << 1, the current becomes linear 
in h and proportional to the (persistent) current that results from a real vector potential. This real persistent current 
was analysed, for arbitrary (7, by Gogolin and Prokofe^i^ and we may simply take over their results, replacing the 
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dimensionless magnetic flux, (p, by ihL. In the case of a relevant pin, g < 1, Eq. (18) of 4(| gives: 

ih 



J cx 



Li/9- 



(4.28) 



We see that the linear response current vanishes as a non-trivial power of whereas it is independent of L for .g = 1. 
As observed in j4^ the L-dependence of the current can be understood from the renormalization group behavior of 
the effective weak tunneling matrix element, F, given by Eq. H4.2()|l . 



r(L) cx L^-i/^'. 



(4.29) 



In the weak tunnelling limit of the non- interacting case {g — 1), the transmission coeffient Tp (x \T\'^ and the persistent 
current J cx \/Tp cx F. For g < I, replacing F by F(L), leads to H4.28|) . In the other case of an irrelevant pin, g > 1, 
the transmission coefficient renormalizes to 1 at low energies and long distances so we expect to recover the result for 
the system with no pin, namely 



J 



ih 



-no- 



(4.30) 



(Note that we require not only 1/h >> L but also that L is much greater than the characteristic length scale required 
to send eo ^ 0. ) 

In fig. 1^1 we show DMRG results for the imaginary current in the case of a relevant pin for system sizes up to 
L ~ 256. The finite-size scaling (|4.28() of the linear-response current is confirmed with high accuracy in the limit 
hL < 1. 




1 

h*L 



FIG. 9: Finite-size scaling of the current (DMRG results) for filling no — 0.25, eo — 2 and a relevant pin {g ~ 0.72). Note the 
data collapse in the linear-response regime hL 0. 



We now consider the other limit, L » ^ = l/2gh, where we can characterize the current in terms of the pinning 
number, Np. In the case, g > 1, where the pin is irrelevant, it is reasonable to calculate the correction to the current 
due to the pin (which can be expressed in terms of the pinning number) in lowest order perturbation theory in eg, 
using the bosonized 6'-representation. We again start with the Lagrangian density of Eq. I|4.1|l . 

In second order perturbation theory in eo, the leading correction to the ground state energy is proportional to: 

/oo 
dT < cos[27r6'(0, t)] cos[27r6l(0, 0)] > . (4.31) 
-C30 

Note that this ground state energy is given by the logarithm of the classical partition function associated with Eq. 
()4.1|l . Upon taking into account the shift of Eq. H4.21|l . we find: 



6Eo cx el 



dT 



cos[2gchi 



(4.32) 



where To is a short distance cut-off. For 2 < 2^ < 3, Eq. (|4.32|l leads to: 



SEq cx el 



-h^s- 



constant • r, 







-23+1 



(4.33) 
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The correction to the current due to the pin, is thus: 



L ah L 



cx ^~eoih)^. (4.34) 



The last entry expresses SJ in terms of the renormahzed pinning strength at the length scale 1/h. After setting g — 1, 
we find agreement with the exact result in the dilute limit, if we also assume that h « ttuq. We only expect the 
phonon representation to be valid when h << tthq since no represents an effective ultra-violet cut off and we have 
thrown away higher derivative terms in deriving this effective Lagrangian. For 2g > 3, the leading behavior at small 
h, from Eq. 1^^, is: 



1j 



Thus the pinning number behaves as: 



« )^ (2 < 2.9 < 3) 



cx {2g > 3), (4.36) 

results applicable in the limit Lh — )■ oo. 

We may straightforwardly extend this calculation to finite Lh. Then we must use the finite L version of the 
correlation function in Eq. (|4.31ll , which can be obtained by a conformal transformation and effectively replaces r by 
(L/ttc) sinh(7rcr/L) in Eq. H4.32(l . Upon rcscaling the r integration variable, u = ttct/L, and using Eq. (|2.35(l . we 
obtain: 

SEo (X el{L/TTc) , - " , — nrC0s[2ghLu/TT], (4.37) 

Differentiating with respect to h and dividing by h gives the pinning number: 



oo 



iVp oc.lL-^^/.-W ,, "^in(2gW.) ^ (l<g<3/2). (4.38) 
Jo smh ^ u 

Note that we have set the lower limit of the integral to since it converges. Thus we obtain the scaling form: 

Np = L^'-^Sf^Lh)^ (4.39) 

where: 

/'^)4r-=is^- ("«) 

We see that f{x) x^9~^ as x ^ 0. Using DMRG, we have numerically checked this scaling for the case 2g < 3 and 
have found good agreement (see fig. ITUI) . 

To study the case of a relevant pin, g < 1, we use the weak tunnelling model of Eq. H4.f 7|l . We determined the 
effective F parameter by using the h = RG equations of Eq. H4.20|l out to a length scale of 1/h: 

Fcx/i^/f-^ (4.41) 

We noted in Sec. Ill that it is convenient to make a similarity transformation of the Hamiltonian such that the 
hopping terms have the form of Eq. (|3.30|) in which all the non-Hermiticity resides on the weak link. We then showed 
that the exponentially small hopping term, oc Ve^^^ from right to left can be dropped at small F, keeping only the 
exponentially large hopping term, oc Te^^ from left to right. Thus we have a "one-way hopping model". It seems 
plausible that this approximation can also be made away from the (? = 1 free fermion case, as long as the effective 
F is very small. Once we make this approximation, F and h appear only in the combination Ve^^ . Then, since we 
expect that, in the infinite L limit the current goes as: i—h, the finite size correction at small F must go as: 



r no 
m 



InF 



(4.42) 
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FIG. 10: Scaling function for the pinning number in the case of an irrelevant pin {U — 10, V = 0, g ~ 1.19) with boson density 
no = 0.25 and pinning potential eo ~ 2. 



We now replace F by its value in Eq. (|4.41l) which leads to: 



i— 
m 



1 In + constant 



(4.43) 



Thus the pinning number is: 



N„ ^ N, 



J(eo = 0) - J(eo) _ no 
J(eo = 0) ~ h 



1 

- - 1 
9 



\\nh\ + C 



(4.44) 



where C is a non-universal constant. Note that the universal number g appears here as an amplitude, not as an 
exponent. We emphasize that we expect this amplitude to be universal by this argument, although the sub-dominant 
constant is not. Thus a relevant pin only increases the pinning number by a factor of |ln/i| compared to the free 
fermion (low vortex density) case where the pin is marginal. 

Our DMRG results for Np in the g < 1 case, shown in fig. ^2 confirm these analytic predictions. The logarithmic 
behavior (|4.44|) is clearly observed for two different values of the Luttinger liquid parameter g. Note that at Lh cx 1 
the divergence of Np is cut off, analogous to the free-fermion limit. 




FIG. 11: DMRG results for the pinning number. Main plot: filling no = 0.25, eo = 2 and a relevant pin {U = 10, V = 4, 
corresponding to jr ~ 0.72). The dashed line gives the logarithmic behavior in Eq. (14.441 with C = 0.5. Inset: same for U = 20, 
y = 10, 3 fa 0.62 and offset C = 0.4. 



In Appendix F we briefly discuss the phase-boson representation of the weak tunnelling model after the similarity 
transformation, Eq. H3.3()|l . for general g, pointing out its connection with a model of current interest in string theory. 
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V. POINT DISORDER 



This paper focuses on the physics of thermal fluctuations in very clean (l+l)-dimensional superconducting slabs. 
The appropriate physical conditions necessary to neglect weak point disorder due to oxygen vacancies, proton irra- 
diation, etc^ in high T^, superconductors are discussed in Appendix E. We can, however, get some insight into the 
influence of a single extended defect in the presence of strong point disorder using the linear response formalism 
discussed earlier. A treatment of a single columnar pin in the presence of strong point disorder, including a full 
renormalization group analysis, will appear in a future publication. 

The effect of a quenched random distribution of point pins on vortex arrays in thin superconducting slabs has been 
discussed in a number of publications^' ^^i^^ . In the presence of point disorder, the background elastic free energy 
H2.14(l on which we impose a single columnar pin becomes 

f 1 1 

—h{dTu) + Vo{x, t) cos [2ttu{x, T)/ao — (3{x, r)] }. (5.1) 

As pointed out in Rcf. 6', the statistical mechanics associated with Eq. (|5.1() is similar to that of an anisotropic two 
dimensional random held XY model with, however, no topological defects such as the dislocations shown in Fig. 
As in charge density wave physics^^, the cosine term acts to locally flx the phase of the (l-l-l)-dimensional vortex 
crystal at P{x, r) with a strength determined by the amplitude Vo{x, r). The couplings ^ix{x, t) and firix, r) are zero 
mean random variables which describe local variations in the preferred density and tilt of the vortex lines induced 
by the particular configuration of point disorder. At long wavelengths, ^x{x,t) and plt{x^t) renormalize in the same 
way and can be taken to be Gaussian random variables described by a single variance c&^, 



lJLi(x,T)nj{x' ,t') = a5ij5{x - x')5{t - t'), (5.2) 

where i,j — x,t and the overbar represents a quenched average over the disorder. We also take Vo(x,r) to be given 
by a Gaussian distribution, with variance 



Vo{x,T)Vo{x',T') = AoSix-x')S{T-T') (5.3) 



while /3(x,t) is uniformly distributed on the interval [0, 27r]. 

When averaged over point pinning, the linear response equation H2.5()|l becomes 

,Vd{x')^ 



Sn{x, t) = - I dx'dr' ^ ^ C(x ~x',t- t'), (5.4) 



with 



C{x — x' ,T ~ t') — < n{x, T)n{x', t') >o — < n{x, t) >o < n{x' , t') >q. (5.5) 
The Fourier space version of these relations reads 

5n{qx,qr) = S{qx,qT)^^^^^^VD{qx), (5.6) 

where 



S[qx,qr)^<Hqx,qr)V>, (5.7) 

and all averages are evaluated in the absence of the columnar defect. 

We first set h = and a pply the renormalization group analysis of the bulk statistical physics problem defined 
above to evaluate S{qx,qT) and hence determine the change in vortex density due to a single columnar pin. The 
recursion relations for the scale-dependent couplings g{l), and Ao{l) arei^i^2i^2i^ 

MO = 0, (5.8) 



(5.9) 
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FIG. 12: Renormalization group flows in g and Ao for the random field XY model which describes (1+1) -dimensional vortex 
arrays subjected to point disorderf— The Luttinger liquid fixed line at Ao = is is stable to point disorder for g > 1, which is 
also the regime where a single columnar pin is irrelevant. This line becomes unstable to a new fixed line where point disorder 
is important for g < 1. 



^ 2[1 - .9(/)] Ao(/) - C2^l{l) (5.10) 

where Ci and C2 are positive constants. The renormalization group flows in the (5, Ao)-plane are shown in Fig. 1121 

The line of fixed points at Aq = in Fig. ^] is stable to point disorder when g > 1. This is also the regime where 
columnar pins are irrelevant in pure systems. However, as the point disorder strength tends to zero, it generates via 
Eq. H5.9|l a nonzero variance for the couplings fixix^r) and ijLt{x^t) in Eq. (|5.1|) . After setting the cosine coupling to 
zero in Eq. IjS.lf) it is straightforward to show^^ that the leading singular term (to = 1) in Eq. H2.48|l is replaced by 

cos(27moa;) 



<n(a;,r)n(0,0)>-n^^p^-5-^, (5.11) 
with 

^-2.9 + ^, (5.12) 

where CToo — hm;^oo "'(O is a (positive) correction to the Luttinger liquid result. Hence, r/ > 2 for all g > 1 and we 
conclude that S{qx, Qt) never diverges near any of the reciprocal lattice vectors. In particular, S{qx, Qt) never diverges 
at Gi, suggesting that an isolated columnar pin becomes even more irrelevant for in the presence of strong point 
disorder when g > 1. 

A new stable line of fixed points appears in Fig. E| for g < I, signaling the onset of a "vortex glass" phase in 
(l+l)-dimensionE^. The analysis of Hwa and Fishei4^ (see also, R.ef. l5]Jl shows that 



(ni(x,T)n*(0,0)) -exp[-C3ln2(a;2+cV2)] (5.13) 



where C3 is a positive constant. Because this decay is faster than any power law, S{qx, qr) is again finite everywhere 
along the g^^-axis, suggesting that an isolated columnar defect is also asymptotically irrelevant in the presence of strong 
point disorder in this regime. However, if the point disorder is weak, the effective pinning strength of the columnar 
defect can grow quite large over intermediate length scales. The subtle and complex physics which distinguishes the 
response of (l+l)-dimensional vortex arrays to a columnar defect above and below the vortex glass transition will be 
discussed in Ref. (j4S) . 

We conclude this brief discussion of point disorder with two comments: The tilt field ft, is a strongly relevant variable 
in (l+l)-dimensions (the recursion relation for h readsi^ dh{l)/dl = h{l)) for both pure and disordered systems. When 
point disorder is present, the structure function S{qx, qr) can only become less singular along the g^^-axis as a result 
of the afSne transfomation discussed above for pure systems. Second, the slowly decaying translational correlations 
which produce interesting physics for vortex arrays in (l+l)-dimensions also appear in the "Bragg glass" phase which 
arises when the Abrikosov flux lattice is subjected to weak point disorder in (2+1) dimensional^. A study of related 
non-Hermitian Luttinger-liquid-like phenomena when a single twin plane or grain boundary is inserted into bulk 
vortex arrays with a tilted field and subject to point disorder is currently in progress>^ 
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VI. CONCLUSIONS 



We have studied interacting vortices in a thin platelet, in the presence of a single columnar pin, using a combination 
of exact methods for the dilute limit based on the free fermion representation, field theory methods and numerical 
techniques. The response to a pin is controlled by the Luttinger liquid parameter, which depends in a complicated 
way on the details of the inter- vortex interactions and on the density. When the pin is parallel to the magnetic field, 
it is an irrelevant perturbation for g > 1 but a relevant perturbation for g < 1. In both cases a single pin produces a 
local vortex lattice with density oscillations that decay as a power law. A transverse magnetic field introduces a new 
length scale, 1/h. The density oscillations decay exponentially beyond this distance from the pin. We characterize 
the imaginary current, or transverse Meissner effect, by a "pinning number", Np, which measures how many vortices 
are stuck in the vicinity of the pin. This number diverges as {l/h) \ \nh\ for g < I but as h?^^^ for g > 1. When g = 1, 
the free fermion mapping shows that Np oc 1/h. Even for zero tilt, point disorder drastically modified the critical 
behavior associated with g passing through 1. However, it may be possible to experimentally probe non-Hermitian 
Luttinger liquid physics using sufficiently clean high- Tc thin platelets with notches cut on the surface. 
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APPENDIX A: CRITICAL EXPONENTS OF THE DILUTE BOSE GAS 



It is well-known that the low energy long distance properties of a Bose gas in the dilute limit (average inter-particle 
separation, I/uq, large compared to interaction range) are given by a gas of free fermions. In the continuum elastic 
theory, defined in Sec. II, this corresponds to the value g — 1 oi the dimensionless Luttinger liquid parameter which 
determines all critical properties. In this appendix we derive a general result for the leading correction to this limiting 
value of 5, expressing our result in terms of the scattering length, a, determined by the boson interaction, V{x) in 
Eq. (|2.2Q(I and the density, uq. Apart from its applications to flux lines in thin platelets, we expect that this formula 
will be of quite general applicability to various one-dimensional quantum models. 

Our result depends critically on two well-known features of the Bose gas. One of them is Eq. (j2.35|l . eg /it — nom. 
This is expected to be exact for a Galilean invariant gas.^"^ This follows by considering the energy of the entire system 
lowest energy of momentum P. The conserved momentum density is; 



ax \ax 



no-f. (Al) 

da; 



Thus the lowest energy state with a very low momentum momentum P is one in which (j){x) — (Px/nQL). From the 
Hamiltonian in (/)-representation, Eq. (|2.34|l . we see that the energy of this state is: 

On the other hand, from Galilean invariance the exact energy of this state, which is simply one in which all N bosons 
are given a boost to a momentum P/N is: 

^oiP) - (A3) 

Comparing these two formulas gives Eq. (|2.35|l . We note that this argument doesn't really require exact Galilean 
invariance; it is enough that the disperson relation be approximately quadratic at small momentum. For instance, 
a quartic term in the dispersion relation would lead to a {d(f)/dx)'^ term in the Hamiltonian in (/)-representation but 
wouldn't interfere with this determination of the coefficient of the {d(f)/dx)'^ term. We also require an expression 
for the compressibility in terms of g and c. To this end, consider the change in energy of the ground state when a 
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relatively small change, SN is made in the number of particles. From Eq. (|2.26() . which expresses a uniform density 
perturbation as Sn — d9/dx and Eq. (|2.40|1 giving the Hamiltonian in ^-representation, we see that 

Therefore the compressibility, k is given by: 

1 _ 7V2 fd^En 



- T \ »AT2 I ^ —< (A5) 
Upon combining Eq. I)A5|I with Eq. H2.35|l . we see that g is completely determined by no, m and n: 



.9 = V^- (A6) 
V m 

We therefore focus on calculating the compressibility k of a dilute Bose gas. 

Some general results on dilute Bose gases in (1+1) dimensions were derived in Ref. (53). There it was argued that, 
in lowest order approximation, the ground state wave-function can be taken to be of the form which is exact for a 
continuum (5-function interaction.^^ If we label the co-ordinates of the N bosons Xi and assume xi < X2 < ■ ■ ■ < xn, 
this result takes the form: 

/ N \ 

'if{xi,X2,--.XN) ~ ^yl(P)Pexp i^kjXj , (A7) 
P V 3=1 J 

for some set of wave- vectors kj ( which must all be different). Here P permutes the fcj's and the sum is over all 
permutations. The A(P) are coefficients to be determined. The many body wave-function, ^E", is determined for other 
orderings of the Xj from the required symmetry of the wave- function which follows from Bose statistics. By considering 
what happens when two particles approach each other, we can see that: 

A(0) = -yl(P)e*2''[('=*-'=^)/2]^ (A8) 

Here the two permutations P and Q differ only by interchanging particles i and j and S{k) is the even channel phase 
shift. The quantity S{k) is defined by the behavior of the 2-particle wave-function at long distances: 

*(a;i - X2) sin[fc|a;i - X2\ + S{k)]. (A9) 

At small k, the limit which concerns us for low density, the phase shift behaves as: 

S{k) -ak, (AlO) 

which we take as the definition of a, the scattering length. The periodic boundary conditions give a set of constraints 
which determine the allowed fc^'s in terms of the phase shift. In the low density limit these conditions become simply: 

kj{L- Na) +a^ks — TTTlj, all j, (All) 

s 

where the integers nj must be all even for N odd and all odd for N even. The solution of these equations can be 
further simplified in the low density limit, N/L << a. In lowest order we obtain simply: 

kjQ = TTHj/L, (A12) 

where the rij are even for N odd but are odd for N even. We may calculate the energy by simply summing the kinetic 
energy in a region of parameter space where the particles are all far apart: 



3 



so we see that we must choose the smallest possible kjo^s to get the ground state corresponding to a "Fermi surface" 
for ID bosons. The Fermi wave-vector is determined in the usual way: 
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By expanding Eq. (|Alip to next order we obtain kj = kjo + Skj, with 



L 

s 

Thus, in the hmit L — > cxd, the ground state energy density to next order is: 

1 r'"' dk,^ a f^" dk f^" dk 



5kj = ^fcjo ^ 7 X! ^"^^ (A15) 



E,/L « — / ^k' + -^ I ^ I ^(fc - k'f. (A16) 



Upon expressing kp in terms of no we obtain: 



Eo/L « + (A17) 
bm 6m 

We can now obtain a formula for the compressibiUty, k, namely 

From K and Eq. ljA6|l we determine the Luttinger liquid parameter g at low density: 

g w 1 - 2ano + ©(a^n^). (A19) 

We note that a quartic correction to the dispersion relation would only correct Eq. ljA17|) and hence 1/k at 0(11^) 
and hence lead to a correction to g of 0(n§). Therefore, we expect Eq. ljA19|l to be correct even for non-Galilean 
invariant systems. Applications of this result to spin chains will be discussed elsewhere^ 



APPENDIX B: DETERMINING g FOR THE TIGHT-BINDING MODELS 

Our predictions about the critical behavior of the system all involve the parameter g. To test these predictions 
via DMRG on the lattice tight-binding model we need to know g as a function of the microscopic parameters of that 
model. At low densities we may obtain an analytic formula for g using the general result Eq. (|l.l(l together with an 
exact formula for the scattering length of our tight binding model. For larger densities we must rely exclusively on 
numerical calculations of quantities, in the tight binding model with no defect and h = 0, for determination of g. 

We now calculate the scattering length for the tight-binding model of Eq. H2.22|l . To find the phase shift we consider 
the sector of the Hilbert space with 2 bosons, total momentum zero and even parity (as required by Bose statistics). 
We may write the eigenstates in terms of the amplitude 'i'j for the 2 bosons to be separated by a displacement j. We 
may also write a lattice Schroedinger equation in this subspace. For \j\ > 2 this equation reads: 

-2<[*,+i+^',_i] =S^',. (Bl) 

The factor of 2 arises because we can increase the separation of the bosons by hopping either one of them; each process 
contributes a term t. The equations are different for |j| = 1 or and read 

- 2i[^'2 + -^o] + V^-i = 

-2t[*i + ^-^i] + [/^-o = -B*o- (B2) 

We write the wave-function as 

^j=sm[k\j\+Si (|j|>l), (B3) 

in agreement with the definition of the phase shift in Eq. (|A9|) . Due to the short-range interaction, this ansatz exactly 
satisfies the Schroedinger equation, with the correct choice of 'J'g and with E = —At cos k. Upon substituting in Eq. 
HB2|I . we obtain: 



- 2i[sin(2A: -f (5) + *o] + ^^sin(A: -f ^) -4i cosfcsin(A: -I- (5) 
-4tsin(/c + (5) + (C/-|-4tcos/c)^'o = 0. 



(B4) 



32 



0.2 



+ 

V 



0.04 



' 1 ' 

- 1 




' 1 ' 1 

1 ' 




-U=10 






— U=20 






-U=30 












.1.1.1.1.1.1 



20 



40 



60 



80 



100 120 



FIG. 13: Bosonic correlation function for the lattice model 12.221 1 with L — 128, V — 0, no — 0.25 and h = 0. Symbols 
represent DMRG results for the bosonic correlation function and lines are fits based on the conformal field theory prediction 
with g = 1.19/1.09/1.06. 



A little algebra gives: 



Upon taking the limit k 



~&t^ sin k + {U + 4:t cos k) [{V + U cos k) sin k ~ 2t sin 2k] 
^ 8^2 cos fc + ([/ + 4t cos k)[-{V + 'it cos k) cos k + 2t cos 2k] ' 

we obtain 6 —ak with: 



AtV - UV 



2tU + UV + AtV 



(B5) 



(B6) 



Note that a diverges as [/, F at fixed t as for the continuum ^-function potential. Thus, at low density, the 
Luttinger liquid parameter, g, behaves as: 



g -^1 + 2no 



8t^ - AtV - UV 
2tU + UV + AtV' 



(B7) 



We emphasize that this result is expected to be valid for arbitrarily large positive U and V at low enough uq. 

We now turn to estimating g from numerical results on finite systems. The most straightforward way of doing this 
is from the correlation function of n or '0 using Eq. 12.48|) or (|2.54|) respectively. Since our numerical results are for 
finite systems with periodic boundary conditions, we instead use the finite size versions of these formulas which can 
be obtained by a conformal transformation: 



< n{x)n{0) > 



constant x 



cos(27rnoa;) 
|Lsin(7rx/L)|25 



constant 
|Lsin(7ra;/L)|i/2s 



(B8) 
(B9) 



Alternatively, g may be extracted from the finite size spectrum in various ways. The compressibility is easily 
obtained from the energy to add or remove one particle: 



1 

K 



CTT 



9 



[Eo{N + 1) + Eo{N - 1) - 2Eo{N)]. 



(BIO) 



However, we need another result to determine separately the phonon velocity c and the Luttinger liquid parameter 
g. A simple possibility is to measure the speed of sound from the excitation energy of the lowest excited state with 
wave-vector k, choosing the smallest possible non-zero wave- vector, k = 2tt/L. 



E{N,k) - Eo{N) w c\k\. 
Another possibility is to measure the 1/L correction to the ground state energy: 

Eq « eoL — — , 



(Bll) 



(B12) 
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FIG. 14: DMRG result for the density correlation function with [/ = 10, V = 4, L = 128, no = 0.25 and h = (black circles). 
The red line is a fit based on the CFT formula IIBSt with g — 0.72. Note the strong Friedel oscillations due to g < 1. 
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FIG. 15: Luttinger-liquid parameter g calculate for the lattice model 12.2211 with ft = eo = as a function of the boson 
density no. We have plotted results for two different values of the interaction parameters corresponding to g > 1 and g < 1, 
respectively. The dashed lines show the analytic result <B7l l at low densities. 

for a non-universal constant, cq. This follows from Refs. I|5Q) using the fact that periodic boundary conditions on ip 
correspond to periodic boundary conditions on (mod 27r) or 9 (mod 1). 

Within our DMRG simulations we have estimated g using the boson correlation function in eq. p9ll . see Fig. 
We also show the fit for the density correlations in Fig. (|14|l . The agreement is very good except at short distances 
where the field theory predictions are expected to fail. We have checked that the values of g so obtained are in very 
good agreement with the values obtained from the compressibility ljB10|l and the 1/L correction to the ground state 
energy (|ljl2|) . Note that g for the Bose-Hubbard model has been calculated previously by Kiihner et al.^^ where 
the authors used periodic boundary conditions in the DMRG algorithm. Our results for g are shown in fig. H15|l . 
For small densities the agreement between the asymptotic expression (jB7|) and the numerical data is obviously very 
good. We expect that with increasing density, g behaves in a non-monotonic fashion and finally approaches .9 = 1 
again as the filling becomes commensurate in the limit no —^ 1. For finite next-neighbor repulsion > and n — 0.5 
a charge-density wave instability occurs^ which we have not studied here. For the main motivation of our work - 
interacting vortex physics - the Bose-Hubbard model is only applicable in the limit of small filling, since the lattice 
constant has no direct physical meaning. Most of the results in our paper have been obtained for no = 0.25. 

APPENDIX C: BOSON CORRELATION FUNCTION WITH A PIN 

In this appendix we discuss numerical and analytic results on the correlation function of the boson creation operator 
across a pin, in the absence of a tilt field, /i = 0. This correlation function clearly distinguishes the cases of a relevant 
((7 < 1) and irrelevant (5 > 1) pin. 
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We note that the finite L equal time boson correlation function at the broken chain (F — > 0) fixed point, for y > 0, 
is given by a conformal transformation of Eq. I|4.13|l : 

/ ,t/ ^,/ / sin(7ra;/£) sin(7ry/£) V^^^ 

1^ sm [t:{x — y) /2L\ sm [t:{x + y) /2LJ J 

This correlation function is now small but non-zero for two points near but on opposite sides of the pin, x « L and 
L — X « L. These points are correlated by going around the chain, rather than across the pin. Even for finite eq we 
expect this result to hold in the case of a relevant pin, 5 < 1, at sufficiently large L, x, y, N — x and L ~ y. We note 
that neither x nor y can be too close to the pin. 

We now consider specifically the boson correlation function across the pin, < ^p^x)^^p{-'x) >, {x > 0) in the limit 
of infinite system size. If the pin is irrelevant, g > 1, then at large x we should recover the result for the system with 
no pin: 

, + . X , / N constant 
< ^\x)^/j{~x) >^ -^TJIT- (^2) 

On the other hand, if the pin is relevant, g < 1, this correlation function vanishes more rapidly with x. The weak 
tunnelling, F, model is again useful for this calculation. We may calculate this correlation function in lowest (first) 
order perturbation theory in F and then "renormalization group improve" the calculation by replacing F by its 
renormalized value at scale x. Perturbation theory gives: 

< V^(x)V(-a;) >^ -r [ dT< V'^(a;,0)V'(0+,r) >o< ^^(0", T)V'(-a;, 0) >o +c.c. (C3) 



The correlation functions inside the integral, evaluated using the Hamiltonian of Eq. H4.17|l with F = 0, may be 
obtained from Eq. (|4.13|) by setting y to a value of order a short distance cut off: 



<^\x,OMO+,t)>o^\j-^^^^^^\ (C4) 



i/4g 

Upon inserting this expression into the integral of Eq. ljC3|) , we obtain 



From Eq. (|4.2()|l we see that 



so finally 



<^\xM-x)>^-^^. (C5) 



reff(^) « -ZTJ^' (C6) 



<^l;\x)i:{-x)>^-^^. (C7) 



Since we are assuming now that g < 1, we see that this correlation function drops off more rapidly than l/x^^"^^ , the 
result for no pin. 

Finally, if we take into account the finite L periodic boundary conditions in the case of a relevant pin, 
< {x)ip{—x) > would be given by a sum of the result for an infinite pin of Eq. IIClll and the result of Eq. 
(IC7|I . suitably generalized to take into account the finite L. For x « L, this is the sum of two small terms. One 
is the weak correlations across the pin and the other is the correlations from going around the circle the long way 
without crossing the pin. 

Numerical results for the boson correlation function are shown in fig. 1161 Evidently, correlations across the pin are 
strongly suppressed if the defect is relevant {g < 1). From the numerics it is unclear whether the exponent of the 
power law decay changes as given by (|C7p . On the other hand, a pin of the same strength has almost no effect if 
g>l. 
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FIG. 16: Impurity effect on the bosonic correlation function. For a relevant pin with U — 10, V = 4 and g = 0.72 (circles) 
we find a strong suppression of correlations close to the impurity (which is at a; = 0). On the other hand, for an irrelevant pin 
with U = 5,V = 0, g = 1.41 (stars) the dependence on the scattering potential eo is very weak. (These are the upper curves 
which lie almost on top of each other.) 



APPENDIX D: PERIODIC/ ANTI-PERIODIC GROUND STATE ENERGY DIFFERENCE 

In this appendix we discuss the difference in ground state energies with anti-periodic and periodic boundary condi- 
tions on ■(/;. In the tight-binding model this corresponds to changing the sign of the hopping term on any single link. 
From Eq. (|2.27|) we see that anti-periodic boundary conditions on -0, 

i;{x + L) = -ij{x), (Dl) 

imply that the phase field (p must wind by a half-integer: 

(l){x + L) = 0(.t) + {2p + 1)7T, (D2) 

for integer p. The ground states correspond to (j){x) = zLirx/L with energy: 

Er < = ^- (D3) 

Upon noting that Eq. IID3|) determines eg while the compressibility gives g/c, we see that both parameters are then 
determined. 

The effects of a pin on the periodic/anti-periodic ground state energy difference can be calculated using "renor- 
malization group improved" perturbation theory. Antiperiodic boundary conditions in the tight-binding model with 
a pin are defined by the Hamiltonian of Eq. H2.22|l but with the periodic condition of eq. H2.23|l replaced by: 

bL = -bo. (D4) 

When the pin is irrelevant, g > 1 and we expect to obtain the result of Eq. (|D3|I asymptotically for large L. On the 
other hand, if the pin is relevant, g < I, the two sides of the system become asymptotically decoupled and this energy 
difference scales to zero more rapidly with L. Note that Eq^ — Eq is strictly zero at the open chain fixed point. 
Within the effective Hamiltonian of Eq. (|4.17() . we see that when F = 0, it makes no difference whether we impose a 
boundary condition ip{0^) — ±-0(0^). For small F, this energy difference is first order in F. From Eq. IjClfl . with x 
and y of order a short distance cut off, we obtain: 

If we now replace F by its renormalized value at scale L, we obtain: 

^r-E^^j^. (D6) 



Since g < 1 here, we see that the decay is always faster than with no pin. 
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In fig. (|17|l we show numerical results for the boundary sensitivity AEq = E^^ — E^ normalized to the value without 
any defect. Although our system sizes were not big enough to extract power law behaviors, the data suggest that for 
a relevant pin {g < 1) the normalized sensitivity decreases to zero in the thermodynamic limit, indicating that the 
system is indeed asymptotically "cut in two" . (For a very large system, we expect A£'o(impurity/A£'o(iio impurity) oc 
L^~^/^.) For g > 1, on the other hand, our numerical results are consistent with a ratio which approaches unity for 
large L, indicating the irrelevance of the pin at large length scales. 
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FIG. 17: Sensitivity of the ground state energy to a change from periodic to antiperiodic boundary conditions for no — 0.25 
and eo — 4. 



The Hamiltonian version of the (^-representation is convenient for studying the periodic/anti-periodic ground state 
energy difference in the presence of a tilt field, h (but no pin). [See Eqs. H2.20|l and H2.45|l .] As discussed above Eq. 
(ID2|I . for anti-periodic boundary conditions, (j) must wind by ± half a period in the ground state: (f>{x) = ±ttx/L. 
Note that the non-Hermitian term does not change the ground state but simply adds an imaginary term to its energy. 
We now have: 

E^P -E^ ^±'-^+7Tcg/2L. (D7) 
The imaginary term is 0(1), unlike the real term which is 0{1/L). 



APPENDIX E: STRENGTH OF POINT DISORDER 



In this appendix we estimate the conditions such that the effect of residual point disorder on the thermally excited 
vortex lines, which are the main subject of this paper, can be ignored. This neglect can be justified only over a limited 
range of length scales. As discussed in Sec. 5, point disorder eventually changes the physics for any value of g. Our 
analysis is mostly based on Refs. (0) and H^), to which we refer readers for additional details. 

The discussion proceeds in two stages. First, we demonstrate that in (l-l-l)-dimensional vortex arrays, typical 
pinning energies per unit length associated with an isolated columnar pin or "notch" greatly exceed the energy 
available from ignoring the pin and simply following an optimal path through the point disorder. Thus the columnar 
pinning energy greatly exceeds the collective effect of point pins. Second, given that at least one vortex line is 
strongly pinned on a columnar defect, we provide a rough estimate of a sufficient condition for thermal fluctuations of 
vortices to dominate over point disorder-induced wandering until one reaches the length scale of the renormalization 
group equations of Sec. 5. Our main conclusion is that (l-l-l)-dimensional "platelet" samples of high-Tc materials 
(disordered only by oxygen vacancies) within 5-10% of the critical temperature could provide a good opportunity to 
study the physics of thermally excited vortex lines discussed in this paper over a range of length scales. 

We first compute the average pinning energy of a single vortex, confined by its neighbors in a (l-fl)-dimensional 
array, due to random point impurities. We define a characteristic "imaginary time" length scale I* at which transverse 
fiuctuations of a vortex line due to point disorder reach the average vortex spacing oq : 

6xil*)^xJ[-X =ao (El) 
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where C = 2/3 in (l+l)"diniensions^ and 

Xc = TV^iAo, (E2) 



lc='^xl. (E3) 

Here Xc is the length scale out to which vortex wandering is described by a thermal random walk, and Ic < I* is the 
corresponding length in the imaginary time direction. Note that ii is the line tension introduced in (|2.1|) and Aq is a 
correlator measuring the strength of point disorder. The parameter Aq is related to the correlator Aq introduced in 
the continuum model of Eq. (|5.3(l by Ao = Aq/hq. The pinning energy for an imaginary-time segment of length I on 
this path is given hy^ 



Up{n^T{-] . (E4) 



Upon solving for I* and using C = 2/3 the pinning energy per length due to point disorder simplifies to: 

Upjl*) ^ 
I* aoT' 



(E5) 



For quantitative estimates, we assume isotropy of the effective mass in the a6-plane and average over the sample 
thickness along the c-axis to obtain the disorder correlator as^ 



3/2 

(E6) 



where Jpi, is the pair breaking critical current, Jcp is the critical current due to point defects and (f>Q — 2 x lO^^Gcm^ 
is the flux quantum. With sample thickness w — 150 nm, penetration depth A — 150 nm, coherence length ^ = 2 nm, 
temperature T = 50 K, and the estimate (valid for high-Tc superconductors^) Jcp/Jpb ~ 0.01, we obtain a pinning 
energy per length 



Upin 



I* 



7 X 10"^^erg/cm. (E7) 



For comparison, the pinning energy per unit length due to a columnar defect is given byi 

Uo/L = ( w 1.2 X 10~^erg/cm (E8) 
\47rA / 

and is thus clearly dominant compared to point disorder by several orders of magnitude. 

A related question concerns whether disorder affects the wandering of a single vortex line in (H-l)-dimensions 
before it interacts with its neighbors. Physically, we require that thermally excited vortices collide and restart their 
random walk several times before pinning due to point disorder becomes strong. Up to constants of order unity, the 
requirement is that the thermal length scale Xc in Eq. I|E2|I above exceed the vortex spacing ao, 

>ao (E9) 



eiA 

Using the same parameters as above, we find 

Xc = 3.5 X 10-^cm (ElO) 

at T = 50^:. Upon assuming that Jcp/Jpb is roughly temperature independent and noting that close to the critical 
temperature eiA oc X'^iT) « X-^{50K)\t\^^^ we have 



^3. 5x10"'^ cm 
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Here 

t={T- Tc)/Tc. (E12) 

For a typical high- Tc compound hke YBCO with T = 80K and a reduced temperature |t| = \T — Tc\/Tc = 0.1 we thus 
obtain 

Xc « 2.7 X IQ-^cm > ao (E13) 

which is larger than a typical vortex spacing for H ^ (f>o/{wXc) ~ 400G. With these estimates in hand, we can use 
renormalization group recursion relations [see Sec. 5 and Ref. lj48l) ] to determine the length scale at which point 
disorder dominates. For clean high Tc samples close to the critical point, this length scale can be very large. 



APPENDIX F: CONNECTION WITH STRING THEORY 



The bosonized form of the weak tunnelling model, Eq. (|4.17(l . generalized to arbitrary can be simplified by 
introducing linear combinations of the (independent) boson fields to the left and right of the pin: 



(/)±(a;) = [(t)ix)±(j){L-x)]/V2. 



(Fl) 



Despite the non-local nature of this transformation, it yields a local Lagrangian since the interaction (tunneling) term 
occurs only at x = 0. This Lagrangian takes the form: 



± 



L/2 



dx < ^ 
1 27r 



dr 



dx 



ihcgV2 I 



dx 



fcos[V20_(O)] 



(F2) 



where f (x F, the weak tunnelling amplitude introduced in Eq. 14.17f) . We see that (/>+ decouples from both the pin 
and the tilt field. 0_ obeys mixed boundary conditions, (j)-{L/2) = 0, d(l)-/dx{0) = 0. 

Alternatively, we may make the similarity transformation of Eq. H3.29|l before rewriting the Lagrangian in terms of 
the phase field, 0. This is equivalent to: 



implying: 



(j){x) (j){x) + ihx, 



(t>- (x) -> {x) -f ih{2x - L)/^/2. 



(F3) 



(F4) 



This eliminates the term proportional to ihdcj)^ / dx but makes the tunnelling term non-Hermitian. Dropping the field 
(^-(_, we obtain the equivalent Lagrangian: 



27r 



L/2 



dx 



dr 



dx 



F 



^/iLgi:V20_(O) _^ g-/iLg-iV2,^_(0) 



In the small F limit, it is permissible to drop the exponentially small term, c>c Fe , leaving the Lagrangian: 



L 



27r 



L/2 



dx 



dr 



+ c 



dx 



Finally, if we make the formal analytic continuation, 0_(x, r) — s- —i(j)^(x,T) we obtain: 



L 



9_ 

'27r 



L/2 



dx 



1 fd^ 

c \ dr 



dx 



+ £g/iLgy20_(O)^ 



(F5) 



(F6) 



(F7) 



This is a strong coupling limit of a model, known as the time-like boundary Liouville theory, which has received 
considerable attention in the string theory literature lately (in the marginal case, g — I), because of its connection 
with "s-branes", i.e. space-like "d-branes". See, for example, Ref. (;59i) . 
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